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Abstract 

(N ; 

We discuss general properties of a system of heavy fermions (including antiparti- 
cles) interacting with rather light bosons. First, we consider one diagram of The 
fermion chemical potential is assumed to be small, [if < T. Already for the low tem- 
CN I perature, T <C min(Tbi.f, m&), the fermion mass shell proves to be partially blurred 

due to multiple fermion rescatterings on virtual bosons, m& is the boson mass, Tbi.f 

■ (<C rrif) is the typical temperature corresponding to a complete blurring of the gap 
between fermion-antifermion continua, mj is the fermion mass. As the result, the 

Qh| ratio of the number of fermion-antifermion pairs to the number provided by the or- 

dinary Boltzmann distribution becomes larger than unity (Rn > 1). For T > m^(T) 
(hot hadron liquid, blurred boson continuum), m|(T) is the effective boson mass, the 
abundance of all particles dramatically increases. Bosons behave as quasi-static im- 
purities, on which heavy fermions undergo multiple rescatterings. The soft thermal 
loop approximation solves the problem. The effective fermion mass mj(T) decreases 
with the temperature increase. For T > Tbi.f fermions are essentially relativistic par- 
ticles. Due to the interaction of the boson with fermion-antifermion pairs, mt(T) 
decreases leading to the possibility of the "hot Bose condensation" for T > T c b- The 
phase transition might be of the second order or of the first order depending on the 
species under consideration. We study in detail properties of the system of spin 1/2 
heavy fermions interacting with substantially lighter scalar neutral bosons (e.g., Na 
system). Correlation effects of higher order diagrams of are evaluated resulting 
in a suppression of vertices for T > m£(T). The abundance of high- lying baryon 
resonances proves to be of the same order, as the nucleon-antinucleon abundance, 
or might be even higher for some species. Further we discuss the system of heavy 
fermions interacting with more light vector bosons (e.g., Nuj and Np) and then, 
with pseudo-scalar bosons (e.g., Nir). For the fermion - vector boson system corre- 
lation effects are incorporated by keeping the Ward identity. In case of the fermion 
- pseudo-scalar boson system correlation effects are rather small. Finally, we allow 
for all interactions. We estimate Rn ~ 1.5 for T ~ Tbi.f proves to be near 

T c b; both values are in the vicinity of the pion mass m n . 
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1 Introduction 

In heavy ion collisions at ultra-high energies a huge amount of secondary 
particles is produced, cf. [1]. At ~ 130 GeV at RHIC the ratio of produced 
negative pions to protons is ir~ /p ~ 9.5, the ratio of negative kaons to negative 
pions is K~/ir~ ~ 0.15, the ratio of antiprotons to protons is p/p ~ 0.71, and 
at ~ 200 GeV at RHIC n~/p ~ 9.3, K'/tt' ~ 0.16, p/p ~ 0.74. The set 
of available data is summarized in tables in [2] . 

Heavy ion collisions at RHIC create a non-equilibrium fireball consisting of 
strongly interacting secondary particles. Due to multiple collisions of particles 
the system rapidly reaches a quasi-equilibrium. Then, the fireball is charac- 
terized by a temperature T(t, r) in the center of mass frame and by a velocity. 
It is commonly believed that at an initial stage a produced fireball is con- 
structed of quarks and gluons. Then during the expansion of the fireball there 
occurs a process of hadronization. The hadronized fireball contains mainly 
pions and also heavier particles like kaons, a (/o(600 MeV)), p(769 MeV), 
u(782 MeV) and other mesons, and nucleons, A(1116 MeV), £(1192 MeV) hy- 
perons, A (1232 MeV) isobars and heavier baryon and antibaryon resonances 
in a smaller amount. At the so called freeze-out stage resonances decouple, 
their momentum distributions freeze and with these distributions one observes 
particles at infinity. Statistical equilibrium ideal resonance gas model, cf. [3] 
and refs therein, and somewhat different models, e.g. [2], that introduces ef- 
fective non-equilibrium occupancy parameters, being applied to the fireball 
break up stage, allow to fit well observed particle ratios. Typical values of the 
non-strange baryon chemical potential used to fit mentioned RHIC data vary 
in the range ^bar — 20-^-40 MeV, and the break up temperature varies in 
the interval T ~ 140 170 MeV. At LHC (y/s ~ 5500 GeV) one expects a 
tiny baryon chemical potential (~ 1 MeV), cf. [3]. This system can already be 
called the hot vacuum. 

Since secondary particle production rates are very large inside the fireball, 
one may expect an importance of effects, which are beyond the scope of the 
usually used simple approximation of elastic collisions of free particles and 
even of a more involved quasiparticle approximation. To describe a system 
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of strongly coupled resonances one needs to develop an approach including 
both particle widths and the dispersion, taking into account particle feed- 
back effects. Although the kinetic description of resonances within the self- 
consistent so called ^-derivable scheme has been constructed (cf. refs [4] and 
refs therein) it looks very complex. Therefore, in order to understand most 
important signatures of processes one needs further simplifications. 

The problem of the behavior of the heated pion-nucleon vacuum (zero total 
baryon charge) was risen by Dyugaev in ref. [5]. A very intuitive considera- 
tion was sketched in an analogy with the description of the electron-phonon 
interaction in doped semiconductors. Even at zero temperature the tail of the 
electron wave function penetrates deeply into the band gap due to multiple 
electron-phonon collisions [6]. Dyugaev conjectured that in the nuclear prob- 
lem nucleons may play the same role as electrons, and pions, as phonons. 
Within the given analogy, in order to construct a qualitative picture of the 
phenomenon, ref. [5] considered nucleons and antinucleons, as non-relativistic 
particles interacting with pions by the non-relativistic pion-nucleon coupling, 
firNipaq- T7T0. Ref. [5] also conjectured existence of the ending temperature 
for the hadron world, above which the system can't be anymore in the hadron 
state due to anomalous production of fermion-antifermion pairs (cf. Hage- 
dorn picture, its recent application to RHIC energies see in [7]). This critical 
temperature (T dec ) was estimated to be < 250 MeV. 

In this paper we focus on the description of the heated quasi-equilibrium 
hadron liquid having a small or even zero total baryon charge (the hadron 
vacuum). More precisely we will exploit that the net baryon density is much 
less than the density of produced baryons (and antibaryons). This allows us 
to neglect particle-hole effects compared to particle-antiparticle effects, that 
essentially simplifies the consideration. The limit is safely fulfilled for fif < T, 
where /if is the fermion chemical potential. We develop a general relativistic 
approach and match it to the non-relativistic one. We use the term "hadron 
liquid" rather than "hadron gas" to stress crucial role of strong coupling. 
Besides we will see that the pion-nucleon-antinucleon attraction is essentially 
enhanced in the p-wave. Thereby a short range correlation naturally comes 
into play. 

To better understand the physics of the phenomenon we first consider prop- 
erties of an idealized hadron liquid consisting of strongly interacting fermions 
of one kind and bosons of one kind. It is assumed that fermions are essentially 
heavier particles than bosons. We argue for the following qualitative picture. 

There are several temperature regimes. The regime T <C m.m{rn%/rn,f, Tbi.f} 
corresponds to a slightly heated hadron liquid (or the hadron vacuum, if the 
fermion chemical potential /// = 0). Tbi.f is the temperature, at which the gap 
between fermion-antifermion continua becomes completely blurred. Typically 
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Tbu(g) is of the order of ~ m n for relevant values of the fermion-boson cou- 
pling constant g, = 140 MeV is the pion mass. Bosons are almost free 
particles in this temperature regime. Fermion distributions begin to deviate 
from Boltzmann distributions due to multiple collisions of each fermion on 
bosons. This deviation increases with the temperature increase. 

If there exists a temperature interval mf/mf < T <C min{Tbi.f, m^,} (a warm 
hadron liquid, partially blurred fermion continuum), then the fermion mass 
shell is already partially blurred due to multiple rescatterings of the fermion 
on bosons. The quasiparticle approximation for fermions fails, if the fermion- 
boson coupling constant g is rather large (e.g., g ~ 10 for a, uj and p meson - 
nucleon (N) interaction). As the result, fermion distributions become essen- 
tially enhanced compared to the ordinary Boltzmann distribution. For realistic 
hadron parameters the regime of a warm hadron liquid can be realized only for 
pions, not for a, uj and p due to their large masses. However the enhancement 
is not too strong for pions, since g ~ 1 (rather than ^> 1) in the latter case. 

With further increase of the temperature, hadron effective masses substantially 
decrease. For mj ^> T > ml(T) (hot hadron liquid, blurred boson continuum), 
where ml(T, g) is the effective boson mass depending on g, heavy rather rapid 
fermions abundantly produce effectively less massive and slower virtual (off 
shell) bosons (boson tadpole diagram) and undergo multiple rescatterings on 
them, as on quasi-static impurities. Due to the width effect (from multiple 
quasielastic rescatterings) the fermion propagator completely looses the for- 
mer quasiparticle pole shape it had in a dilute medium. For rrif ^> T > T bl f 
fermions become essentially relativistic particles. The hot hadron liquid comes 
to the regime of the blurred fermion continuum. The fermion sub-system rep- 
resents then a rather dense packing of fermion- antifermion pairs. Since also 
T > m* b {T), this state is the state of the blurred hadron continuum (blurred 
continua for both boson and fermion sub-systems). The fermion- antifermion 
density, pjj grows exponentially with the temperature in a wide temperature 
interval. Bosons rescatter on fermion- antifermion pairs (fermion-antifermion 
loop diagram) and due to that decrease their effective masses. At a temper- 
ature T > T cb , the effective scalar boson mass may vanish and a hot Bose 
condensation (HBC) may set in by the second order phase transition. We call 
it HBC, since the condensate appears for the temperature larger than a crit- 
ical temperature. For vector and pseudo-scalar bosons (if the latter interact 
with fermions via pseudo-vector coupling) the HBC may arise by the first or- 
der phase transition at finite value of m* b {T ch ). Moreover, for scalar and vector 
bosons HBC occurs in the s-wave state, whereas for pseudo-scalar bosons with 
the pseudo- vector coupling to fermions the HBC may arise in the p-wave state. 
For realistic values of hadron parameters, the problem of the determination 
of Tbi.f, T c b and ml(T) is the coupled-channel problem. As the result of its 
solution, T c b proves to be close to Tbi.f • In spite of large values of bare masses 
of a, u, p mesons and nucleons, numerically both values T w .f and T cb prove to 
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be in the vicinity of the pion mass m n . At such a temperature the resulting 
density of fermion-antifermion pairs is estimated as < po, where po — 0.5m^ 
is the nuclear saturation density, h = c = 1. 

With subsequent increase of the temperature, for T > T bl f the baryon-antibaryon 
density continues to increase, then it reaches a maximum (the value p max may 
several times exceed po) and then may even begin to drop down. 

The strange particle production, as well as the production of other baryon 
resonances, are significantly enhanced with the increase of the temperature. It 
is known that at high energies of heavy ion collisions the experimentally ob- 
served kaon to pion ratio becomes energy independent. It is usually associated 
with the quark deconfinement, cf. [8,3]. However such a behavior can be also 
naturally explained within the pure hadron picture. In this sense one may 
speak about a quark-hadron duality: observable can be explained in terms 
of both the quark-gluon degrees of freedom, cf. [9], and only hadron degrees 
of freedom [10]. Taking into account multiple scattering effects, the number 
of hadron degrees of freedom is significantly enhanced simulating the same 
effects, as from deconfined quarks. 

In reality hadron and quark-gluon degrees of freedom may interact. Incorpo- 
rating the quark structure of hadrons and the possibility of hadron states in 
the quark matter, one may expect strong thermal fluctuation effects of the 
quark-gluon origin in the hadron phase, and strong fluctuation effects of the 
hadron origin in the quark-qluon state. 1 Thereby, instead of a sharp first or 
second order hadron-quark phase transition one may expect the existence of 
a broad region of a hadron-quark continuity, cf. arguments for the crossover 
from lattice simulations [12]. Thus, more likely, at such conditions the sys- 
tem state represents a strongly correlated boiled hadron- quark-gluon porridge 
(HQGP) rather than the pure hadron or pure quark-gluon state. The pure 
quark-gluon phase probably occurs at a higher temperature, T > (1.5 -j- 2)m n , 
and a pure hadron phase, at T < (0.7 -j- l)m n . Below, discussing a high tem- 
perature regime we artificially disregard the quark-gluon effects postponing 
their study to the future work. 

In the low temperature regime, T <C min(m(,, T^u), fermion energies of our 
interest are p — m/ — — m b — 0{mfT /m b , ^JrrifT), virtual boson energies 

are go — m b + 0(mfT/mi„ ^JrrifT), and fermion and virtual boson momenta 

are of the same order |p| ~ \q\ ~ ^j2mjT. In the high temperature regime, 
T > m* b (T), the quantity J, which we further call the intensity of the mul- 
tiple scattering, is much larger than the temperature squared, if the cou- 

1 This might be in an analogy to that was recently found for the phase transition 
to the di-quark condensate state [11]. The fluctuation region might be very broad 
there. 
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pling constant g is rather high (e.g., g ~ 10 for the Na interaction). Then, 
typical departures of fermion energies from the mass shell are substantially 
higher than those for bosons. Typical fermion momenta increase with the 
temperature from non-relativistic values ~ ^j2Trrif) to relativistic ones 
~ m/). These values are significantly higher than typical boson momenta 

\q\ ~ {^2Tm* b (T),T}. It allows one to consider fermions, as hard particles, 
and bosons, as soft ones, that greatly simplifies the analysis and, actually, 
solves the problem. We call such an approximation the soft thermal loop (STL) 
approximation. 

We argue for a huge stopping power in the course of highly energetic heavy 
ion collisions due to mentioned above multiple collisions. After a short non- 
equilibrium stage the system continues to live rather long at a quasi-equilibrium 
in the center of mass frame undergoing a slow expansion into vacuum. Due to a 
large hadron density and a softening of the vector meson spectrum the dilepton 
production is expected to be enhanced, that can be measured. Distributions 
of particles radiating at the break up stage of the fireball are also enhanced 
compared to Boltzmann distributions at given temperature, since at least a 
part of a large number of virtual particle degrees of freedom concentrated in 
the fireball before its break up can be transformed into distributions of parti- 
cles measured at infinity (however the value of the enhancement depends on 
the dynamical mechanism of the break up, cf. [13,14]). 

The paper is organized as follows. In sect. 2 we introduce a $ approach of 
Baym [15] and its application to the description of the system of coupled 
fermions and bosons. Technical details of the formalism: relations between 
Green functions and self-energies and the scheme of the particle-antiparticle 
separation are deferred to the Appendix A. Sects 3 and 4 are key sections of our 
study. In sect. 3 we construct the description of the coupled fermion - boson 
system within the simplest diagram of ^-functional. The main approximation 
we use is that fermions are assumed to be effectively significantly heavier than 
relevant bosons. We first focus our discussion on the low temperature limit 
T <C (m^jTbi.f) and then on a high temperature limit rrif ^> T > m^(T). 
In sect. 4 we study an example of a fermion - boson system coupled by the 
Yukawa (spin ^-fermion - scalar neutral boson - fermion) interaction. We 
include higher baryon resonances into consideration and evaluate correlation 
effects. The scheme is then applied to other interactions, between fermions 
and vector bosons (sect. 5), and fermions and pseudo-scalar bosons coupled 
by the pseudo- vector coupling (sect. 6). In sect. 7 we schematically discuss 
the behavior of the state of the hadron porridge, when all interactions between 
different particle species are included. In Appendix B we present relevant non- 
relativistic limit formulas for heavy fermions. 

Throughout the paper we use units h — c—1 and the temperature is measured 
in energetic units. 
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2 ^-derivable Approach 



2.1 Dyson equation 

Let us consider a system of interacting fermions and bosons. It is described by 
a coupled channel system of Dyson equations for the fermion (/) and boson 
(b) Green functions 



tions, and Ej are fermion and boson self-energies. All quantities are operators 
in the spin space. For the description of an arbitrary non-equilibrium system 
all the values in (1) are expressed in terms of the non-equilibrium diagram 
technique. For that aim one may use contour, or matrix {— , +} notations, cf. 
[16,4]. For the further convenience we prefer the latter. For the sake of brevity 
we will often suppress tensor indices and sometimes sign {— , +} indices using 
symbolic hat-operator notation. 

2.2 $ functional and self- energies 

As it is known, perturbation theory fails to describe collective phenomena. On 
the other hand, coupled Dyson equations (1) can't be solved exactly and ap- 
proximation schemes are required. Among different approximation approaches 
the self-consistent ^-derivable method seems to be promising. It keeps exact 
conservation laws and exact sum rules. For the quark-gluon plasma it provides 
the quantitative prediction of thermodynamic characteristics, which match lat- 
tice results down to 3 times deconflnement critical temperature, T ~ 3T dcc , 



Assume that fermions are coupled to bosons with the help of a two-fermion - 
one-boson interaction. Then, the "^-functional is given by the series of diagrams 



where the bold solid line corresponds to the fermion/ ant if ermion full Green 
function and the bold wavy line, to the boson/antiboson full Green function, 
small dots denote free vertices. 




cf. [17]. 



i <3> = - 




+ . . . 



(2) 
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Fermion and boson self-energies are obtained by the variation of the functional 
$ over the corresponding Green function (the cut of the line in (2)) 



2 for neutral bosons, 

-^/, 6 = T— tttM , ( 3 ) 

oyiLrffi) i £ Qr f erm i ons anc [ charged bosons. 
The fermion self-energy is determined by the diagram 

^ (4, 

and the boson self-energy is given by the diagram 




(5) 

The fat dot in (4), (5) symbolizes the full vertex. In case of the ^-derivable 
approximation, that deals with finite number of diagrams in (2), the fat vertex 
symbolizes the full vertex related to the given $. 



2.3 The simplest ^-diagram 



Let us restrict ourselves by the consideration of the simplest <£> (the first dia- 
gram (2)). Then the fermion self-energy takes the form 



(6) 



and the boson self-energy reads 




(7) 



All the multi-particle rescattering processes 



(8) 



are then included, whereas processes with the crossing of boson lines (corre- 
lation effects) like 
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(9) 



are not incorporated. 

Using diagrammatic rules and relations of Appendix A between "{— , +}" and 
retarded ("R") and advanced ("A") Green functions (A.l), for the self-energy 
we find 



+ J V (q)Gf(p + q )Vo(-q)G^(q) A 



- / V Q {q)iG- f + (p + q)V (-q)2ImGK(q)^L , (10) 



(2tt)' 

(2tt)* ' 



where we took into account that / V (q)Gf(p + q)V (—q)G^(q)-^^ = 0, V 
is the bare vertex. E.g., for the coupling of the s/ = 1/2 fermion with the 
scalar (sb = 0) neutral boson the interaction term in the Lagrangian density 
is L int = -g s ^(p^ and V = g s . 

Then, separating boson particle and antiparticle contributions (using (+) and 
(— ) sub-scripts in notation of Appendix A) we obtain 



i£*(p) = / Q{qo)%{q)G] + {p + q)%{-q)Gl +) {q)^ 

+ / Q(go)%(-q)^ + (p-q)%(q)Gt ( -)(q)^ 
+ J e(q )V (q)Gf(p + q)V (-q)G-+ +) (q) A 



(2vr) 4 

Ri 



d 4 q 



+ J Q(qo)Vo(-q)Gf(p-q)V (q)G^_ ) (q)j^- 4 
+ / e(qo)Vo(q)tGj + (p + q)V o (-q)A bM (q)-0^ 



-/ e(ft)V (-?)^7 + (p-?)Vo(?)4(-)(9)(2^. (ii) 

O(x) is the step function, and are spectral functions for i = b (or 

i — f) particle and antiparticle respectively, see Appendix A. 
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3 Heavy fermions and less massive bosons within simplest $ 



3.1 Fermion and boson self-energies 



We will continue to deal with the simplest $ (the first diagram (2)). Assume 
that boson occupations are essentially higher than fermion ones. This condition 
is obviously fulfilled for low temperatures T <C m&, since in our case to/ > to&. 
Thus the ratio of fermion to boson occupancies is r oc exp[— (m,f— m^/T] <C 1. 
For a higher temperature T > to|(T), for m*f(T) ^> to|(T"), one estimates r oc 
exp[— (mj — to|)/T] <C 1. For a still higher temperature, mj ^> T > TOj(T), 
there remains a power law suppression, r oc ml(T)/T < 1. Thereby, in the 
whole temperature interval of our interest we may retain in (11) only terms 
proportional to boson occupations. Then eq. (11) yields 



,3 J 

£ ^ Swpl ^^(?)^( P + g )y (-?)A( + )(?)^( + )(?o) 





oo 



+ / I ~ ^)4(-)(?K(-)(<?o). (12) 

From this expression we recover the fermion width (see eq. (A.5)): 





oo 



+ / 7^7^^(-g)^( P -g)l>o(g)4(-)(^)^(-)(*). (13) 
^ ' o 

In eqs (12), (13) in order not to complicate the consideration we omitted 
contributions of quantum fluctuations, which do not depend on boson occu- 
pations. These pieces need a proper renormalization. Within the $ derivable 
approach the renormalization procedure was derived in [18]. In principle, men- 
tioned terms are responsible for important effects. E.g., widths of the A isobar 
and the p meson in vacuum are due to such effects. Here we are interested 
in temperature effects. Therefore we for the sake of simplicity consider only 
thermal contributions assuming that the latter being dominating terms. 

The boson self-energy related to the first diagram of $ (2) is as follows: 



ReEf (q) ~ -2Tr J ^Mq) [ReGf(p + q) + ReGf(p - q) 
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xV (-q)lmGf(p)n f (p Q ), 



(14) 



74 

ff(g)^4Tr | Pv (q)1mGf(p + q)V (-q)lmGf (p) 

x bv(Pt>) -n/(po + 5o)]- (15) 

Weusedeqs (A.9) and that ri/(po+<7o)(l-™/(Po)) = ™&(?o)[™/(Po)-™/(Po+<7o)]- 

We stress that fermion and boson Green functions entering above self-energies 
are full Green functions, although calculated within only one diagram o/$. 

We will further study two temperature regimes of low and high temperature. 
Below we demonstrate that for rrif ^> T > ml(T) a high temperature regime is 
realized. The boson continuum is then blurred. The intensity of multiple scat- 
tering is not exponentially small anymore. For T > Tbi.f the fermion effective 
mass essentially decreases with the increase of the temperature and the den- 
sity of fermion- ant if ermion pairs is enhanced accordingly. Thus the fermion 
continuum becomes blurred. As the reaction on the increase of the number of 
fermion-antifermion pairs for T > Tbi.f, the boson effective mass being signifi- 
cantly reduced compared to the bare boson mass, if initially were m b ^> T b i.f. 
As a rough estimate we obtain TO£(T b i.f) ~ T b i.f for T ~ T b i.f. 

For T <C Tbi.f, the fermion-antifermion density is exponentially small and, 
thus, ml(T) ~ m b . Since ml(T hU ) ~ T bL f, if m b > T w .f, then we also have T <C 
m b . If m b < Tbi.f, the fermion-antifermion density is exponentially suppressed 
for T <C m b . We call the temperature regime T <C min(m 6 , T bLf ) the low 
temperature regime. 

3.2 Low temperature limit (slightly heated and then, warm hadron liquid). 
Urbach law vs. Boltzmann law 

For temperatures T <C min(m&, Tbi.f), for m b , being essentially smaller than 
to/, fermion occupations (oc e~ mf/ ' r ) are much less than boson ones (oc e~ m& / T ), 
and we may use simplified eqs (12), (13). 

In the quasiparticle approximation the fermion spectral function (A. 5) reads 

Af{p) ~ Aj(p) • 2n5 (pi -e 2 p - jTr[A°(p)ReE*(p)]) , (16) 

e p = +^Jm 2 j + p 2 . In the very same approximation the boson spectral function 
becomes 
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* MQ) ■ ^ [q 2 - ul - — Tr[A^(g)Re£^(g)] j , (17) 



q > 0, uj q — ^Jml + q 2 , A°, are spin operators. Index j counts scalar 
bosons (j = s) and vector transversal (j = t) and longitudinal (j = I) bosons. 
A| = 1 for scalar and pseudo-scalar bosons, whereas for vector bosons one 
deals separately with the transversal spectral function (A£ = T^) and the 
longitudinal spectral function (Aj, = L^ u ), see eqs (133), (134) below. A° is 
introduced in (A. 2). 

The quasiparticle approximation is valid, if the particle width is much smaller 
than all other relevant quantities in the given energy momentum region. In the 
quasiparticle term (16) we may neglect the contribution of the ReSj oc e~ m *>/ T . 
In (17) we may omit the fermion particle- antiparticle loop term Re£& oc 
e -2m f /T ( s j nce j 1 <^ mb < mf). After that these spectral functions are re- 
duced to free ones. 

Outside the validity of the quasiparticle approximation, in case of still rather 
small fermion width, from (A. 5) one obtains a regular contribution 



A)^{p) ~ f/(p) (ReiGfip))- 1 )- 2 + 0[Y% (18) 

is the fermion width operator introduced by eq. (A. 5). If integrated, the 
first term should be understood in the sense of the principal value. The con- 
tribution oc 0[T^] gets an additional exponentially small factor and can be 
dropped. Also ReGf (p) ~ ReG°f R (p) in the approximation we use. The full 
fermion spectral function is the sum of the quasiparticle and regular terms 
related to different energy-momentum regions, cf. [19]. The boson spectral 
function still can be considered within the quasiparticle approximation, since 
|ReS fe | |ReSj| for energies and momenta relevant for the low temperature 
case, that we now discuss. Also the consideration is essentially simplified, if one 
describes fermions within the non-relativistic approximation, see Appendix B. 
This approximation is, indeed, fulfilled in the low temperature limit, since 
particle energies are near the mass-shell and typical thermal momenta are 
small. 

First approximation to find the fermion width is to use quasiparticle spectral 
functions (16), (17) suppressing there small self-energy dependent terms, i.e. 
reducing these spectral functions to spectral functions of free particles. Then, 
replacing these spectral functions into (13) we obtain 



f f ~ J e(q )5(q - 



LO. 



dq d 3 q 
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X 



+n b ,(-)(qo) 



<*(P0 - % - Zp-q) 



2e 



V (-q)A° f (p-q)V (-q) 



p-q 



(19) 



Integrating (19) in q we find 



ff(p) 



d 3 q exp[-(u q + Hb)/T] 
(27r) 2 Uqtp+q 



S(p + uj q - e f+q -)S°(q) 



d 3 q exp[-(v q - fi b )/T\ 

(27T) 2 



^qtp—q 



6(po - u q - €f-?)S°(-q), 



(20) 



where S°(q) = V (q)A° f (p + q)V (-q)A%(q) and we used that Ag(g) = Ag(-g) 
and that T <C m&. The contribution of the second term in (20) to the fermion 
occupations is e~ 2mb ^ T times smaller than that of the first term for typical 
fermion energies and momenta. Thus the second term can be omitted. The 
first term corresponds to the energy p ~ m,f — nib < rrif, for T — > 0. 

We may present the fermion 3-momentum distribution (A. 16) as the sum of 
two contributions 



n f ,(±)(p) = n°f, { ±)(p) + 6n ft(±) (p). 
First term is obtained, if one substitutes (16) into (A. 16): 



(21) 



n 



/,(±) 



(p) 



7oA° f (p) 



2e, 



™Bol,(±)(p), 



(22) 



where 71bo1,(±)(p) — exp[— (e p ± jif)/T] are Boltzmann occupations of parti- 
cles and antiparticles. We dropped exponentially suppressed contribution of 
Re£/ oc e~ mb l T . Replacing the quasiparticle width term (20) into eq. (18) 
we evaluate the regular contribution to the spectral function. With the help 
of eq. (A. 16) we find the term 6rif t (±)(p). Since \p\, \q\ <C m,f, we may put 
ep +q - ~ mj, if the latter term does not enter the exponent. In the exponent we 
use the expansion ep +q - ~ e p + pqj e p + <f 2 /(2e p ). Moreover, we will exploit that 
pi — e 2 ~ 2e p (p — € p ). Taking off the integral in \q\ — ^Ju q — mf, using the 5 
function we arrive at the expression 



£p 1 

Sn f ,(±)(p) ^ -^-^ J ioS°[A° f (p)] 2 W*\dPo J dx 

p e v -m h -1 
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xcxp 



T e n T 2e„T 



^2 

(23) 



Here the boson momentum |<f*| ~ J (e p — po) 2 — m 2 . We see that typical mo- 



menta of our interest are \p\ ~ \q\ ~ ymjTCmj. Then integrating the rest 
over the angle we obtain 



%,(±)(P) - £2^^(^)^Bol,(±)(^), (24) 



with 



?.(*) = A /t^IW^^^'V?^. (25) 
4m^ J 1 yr\ z A 

= m b Vz 2 - 1/JepT, i] = \p\/Je^T, z = (e p - p )/m b . 



y 

One can see that the integral is cut off at z 2 ~ 1 + m/T/m 2 (corresponding 
to y ~ 1). 

Eg. f££j i/ie A;e?/ expression of this sub-section. Supposing for a rough esti- 
mate that S° does not depend on q we obtain an estimation 



Sn u±) (p) ~ (//(24 7 r 2 ))(m / T/m fe 2 ) 3 / 2 n Boli(±) (26) 

for T <C ml/rrif, where g is the typical value of the coupling constant. Notice 
that we continue to consider the case T <C min(m fe , T b i.f). We see that the cor- 
rection to the ordinary Boltzmann distribution is small, Snf^±)(p) <C 71b i,(±) 
for T <C m 2 (247r 2 ) 2/ ' 3 /[5f 4//3 m/]. We may call such a temperature regime "a 
slightly heated hadron liquid". 

For min(mb, Tbi.f) ^> T > m 2 /m/ (we may call this regime "a warm hadron 
liquid") we get 



Sn u±) (p) ~ (//(167r 2 ))n B oi,( ± )ln^, (27) 

m b 

i.e., for g ~ 10 of our interest the fermion distribution could be up to sev- 
eral times (depending on g) enhanced compared to the ordinary Boltzmann 
distribution. 

To find the regular contribution to the fermion width one replaces (18) into 
(13). Then we obtain an integral equation 
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fr-/e(,oK(, )%o-^)^g 

X Vo(q)r} e Z(p + q)[ReG°/ R (p + q)] 2 V (-q)A° b (q) 
+ /e(g K(g )%o-^)^^ 

xV (-q)r r f es (p - q)[ReG°/ R (p - q)} 2 V (q)A° b (-q). (28) 



Eq. (28) is much more involved than the quasiparticle term, eq. (19). In the 
perturbative regime (<5?i/,(±)(p) <C riBoi) eq. (28) can be solved iteratively 
yielding small corrections to the quasiparticle estimate. However in the limit 
T > ml/rrif the perturbative consideration is valid only, if the expansion pa- 
rameter, being ~ g 2 /(4uir 2 ), is much smaller than unity. The value v depends 
on the choice of the fermion-boson coupling. Typically, one has v ~ 1. For re- 
alistic values of the meson-nucleon coupling constant (e.g., (7 ~ 10 for a, u, p) 
multiple rescatterings of the heavy fermion on light bosons should be taken into 
account in all orders (for T > ml/rrif ). In the physics of solids a similar effect 
(Urbah law) is well known for the case of the electron-phonon interaction in 
semiconductors. A long tail of the electron wave function arises inside the band 
gap of the semiconductor [6] . For massless phonons the effect is stronger than 
for massive bosons (if couplings in both cases are of the same order of mag- 
nitude). The limiting case T ^> m 2 /mf estimation (a warm hadron liquid), if 
done with an appropriate vertex V , is relevant for phonons. The description 
of a slightly heated hadron liquid is different from that for massless phonons. 

Concluding this sub-section, we have shown that even at low temperatures 
(for "a warm hadron liquid") fermion particle- antiparticle densities might be 
essentially enhanced compared to quasiparticle ones (Boltzmann law) due to 
multiple rescatterings of fermions in the thermal bath of bosons. 

To calculate particle distributions explicitly we need to know the explicit form 
of the spin structure operator 7o5 <0 (g)[Aj(p)] 2 , that depends on the choice of 
the coupling between particle species. Since the fermion energy departs only 
little from the mass shell, \5po\ <C po, and the fermion momentum is also 
small \p\ ~ ^2mfT C mj, from the very beginning one could use the non- 
relativistic approximation for heavy fermions. We did not do it in order not 
to spoil our general relativistic approach, which we further use to describe the 
high temperature regime. 



3.3 High temperature limit ("hot hadron liquid" ). Multiple rescatterings 



Now we will consider a high temperature regime, T > m* b (T). In this case 
the boson continuum is blurred. The temperature exceeds the effective gap 
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between particle-antiparticle continua and the corresponding antiparticle den- 
sity is rather high thereby. As we argue below, in a wide temperature range 
the departure of the fermion energy from the mass shell Spo ~ \fl (see eq. 
(70) below) is much larger than that for bosons, 5q ~ max{m fe — m^(T),T}, 



and typical fermion momenta \p\ ~ J2mf(T)T are much higher than typical 



boson momenta \q\ ~ max{y 2ml(T)T, T}. Therefore we are able to drop a 
g-dependence of fermion Green functions in (12). Then eq. (12) is simplified 



as 



£?(p)~J.G?(p) 

= 1^1^ [^o(?)^?(p)Vo(-?)^( + )(?)n6,( + )(9b) 



(29) 



+Vo(-q)Gf( P )V (q)A bi{ _ ) (q)n b , { - ) (q ) 



where 



JW I 2tt r (9)^/(P)^ (~9)A(+)(9) n 6,(+)(®)) 



(2tt) 3 7 2tt 

v y o 



+H(-g)A / (p)yo(g)A(-)(g)n 6l (-)(g ) A7 1 (p). 



(30) 



If we for a moment ignored a complicated spin-structure, we could associate 
the quantity J with a tadpole diagram 




(31) 



This diagram describes fluctuations of virtual (off-mass shell) bosons. For 
bosons, which number is not conserved, we have fib = and n^+^go) = 
n b,(~)(lo)- We formally presented Gf(p) = Af(p)Gf(p), where the spin struc- 
ture term A/ is separated from the term Gf(p) related to dynamical degrees 
of freedom, cf. eqs (34), (35) below. Then the dynamical part of the fermion 
Green function Gf(p) decouples from the integral. For the case of a hard exter- 
nal fermion having a large 3-momentum compared with the typical momentum 
transfer in the loop we may use the soft thermal loop (STL) approximation. 
The latter is opposite to the hard thermal loop approximation of the soft ex- 
ternal particle with a small momentum compared with the typical momentum 
transfer in the loop. The hard thermal loop approximation is widely used in 
the description of the quark-gluon plasma, cf. [17] and refs therein. In the STL 
approximation we drop the dependence of the internal fermion Green function 
in the loop on the internal momentum transfer. Please notice that only in the 
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high temperature limit T > ml(T) it is possible to drop the g-dependence of 
the fermion Green function. Considering the low temperature limit we were 
forced to retain the g-dependence of the fermion Green function in the cal- 
culation of the fermion self-energy, since there \q\ ~ \p\ for typical values of 
momenta, as it followed from eqs (20), (25). Besides, in the low temperature 
limit in case of a slightly heated hadron liquid we used an expansion of the 
full fermion Green function near its non-perturbed value Gf . In the high 
temperature limit the full fermion Green function is obtained straight from 
the Dyson equation (1). The latter equation is greatly simplified in the STL 
approximation and reads 

[G^ipT'Gf (p) = 1 + J-Gf(p) ■ Gf(p). (32) 

This is the key equation of this sub-section. A perturbative analysis of eq. (32) 
is possible only for J <^ T 2 (see the corresponding estimate after eq. (60) 
below). However the latter limit is not realized within the high temperature 
regime for the case of a strong coupling (g ~ 10), see eq. (66) below. 

The operator A f and the quantity G R are complicated functions of invariants. 
To determine them we present the fermion self-energy in the most general 
form as 

£ / = E 1 ^ + E 2 m / + E 3 ^ (33) 
is the 4-velocity of the frame. The Green function takes the form 

G f = G 11 p / + G 12 m f + G 13 i = A f G f . (34) 
In the rest frame = (1, 0) and (34) is simplified as 

Gf(p) = G 1 (po,p)p , + G 2 (po,P)rn f + G 3 (po,p)pj. (35) 

This equation shows that in general in the rest frame the Green function 
depends separately on p and p. 

The term G 3 yields only small corrections to thermodynamical quantities in a 
wide temperature region. Moreover, the fermion- ant ifermion density, which we 
are interested in below, does not depend at all on the term cx fry in the fermion 
Green function. Thus G3 can be omitted. Also in many cases one may put 
G\ ~ G2 with appropriate accuracy, that further simplifies the consideration, 
since then Gf(p) ~ GiAj(p) is determined by only one quantity (Gi). 

Let us consider spin-zero bosons, Vo = g s . Then the Dyson equation (32) has 
clear diagrammatic interpretation. It can be presented as follows 
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(36) 



describing the fermion propagation in an external field 

■J/gl 



/>ext 



j 



9l 



d 3 q 
(2^f 



dqo 
2tt 



A b ,(+)(q)n b , {+) (qo) + A 6)( _)(g)n 6i (_)(go' 



(37) 



compare with (30). As follows from eq. (29), in difference with the standard 
Dyson equation in the external field the r.h.s. of eq. (36) contains two full 
fermion Green functions and two lines of the external field. Eq. (36) shows that 
the propagating heavy fermion undergoes multiple quasi-elastic rescatterings 
on pairs of quasi-static boson impurities. Impurities are quasi-static in the 
sense of the above used STL approximation. In the given approximation the 
quantity J does not depend on the external frequency and 3-momentum. The 
value J is proportional to the density of impurities. Thus, it demonstrates the 
intensity of the multiple elastic scattering. To better understand this, one may 
compare the non-relativistic spin-averaged limit expression 



7 ^— — — n.rcl 1 

P?(p)] = TTr 



/v 1 2m/ 



J n - rel • [Gf ( P T re \ 



(38) 



cf. the first line of eq. (29), with the quasiclassical non-relativistic equation 
[14,20], 



S7 n - rc1 ^ 4^(20-^(0). 



(39) 



Here is the particle (heavy fermion in our case) non-relativistic self- 

energy, 



-jn.rcl = 1^ Vsty + m f )V s {j) + m f ) j 
4 Amj 



(40) 



Fj n rel (0) is the non-relativistic fermion forward scattering amplitude in the 
medium of independent static scattering centers, 



(41) 
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Pb is the density of centers (in our case the density of quasi-static boson im- 
purities, which we introduce as pb = 2m^|0° xt | 2 ). In more detail different 
non-relativistic limit expressions for fermions are discussed in Appendix B. 

In general case, e.g. for vector bosons, eq. (36) has only a symbolic meaning. 
It is the operator equation for several values of intensities of multiple quasi- 
elastic scattering and for coupled functions Gf , G 2 and G%, which determine 
the fermion Green function. 

The STL approximation may allow to develop a simplified kinetic description 
of the non-equilibrium system with the help of the 3-momentum fermion dis- 
tribution function. Such a kinetic scheme could be then spread out to describe 
coherent di-lepton radiation processes in an analogy to the kinetic description 
of the Landau- Pomeranchuk-Migdal effect, cf. [21]. However these problems 
are beyond the scope of this paper. 



4 System of heavy spin 1/2 fermions and less massive scalar neutral 
bosons. 



To avoid complications with spin-isospin degrees of freedom, as the simplest 
example, we will consider a system of spin 1/2 (sj = 1/2) fermions and spin 
zero (sfe = 0) neutral bosons coupled by the Yukawa interaction, 



Lm = -g s ^4>^- (42) 

In this case A° f = (p 1 + m f ), A° b = 1, % = g s , (G°' R ) _1 = p 1 - m f + iQ, and we 
assume, as before, mj > uib- 

Results of this section can be applied for the description of the Na sub-system. 
In subsequent sections we consider Nu, Np and Nn systems and summarize 
results. 



4-1 Low temperature limit 



From the first term of eq. (20) we find 



f q f p (po - e p ) 



& 2 (7o + i: 

47T 



exp 



Po)/T] 



x^(P0-e p y-m b \ y 



sinh (yy) c _ y 2 /2 

VV 

yj (Po - e P ) 2 - ml \p\ 



(43) 



Po < e p . 
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Using that typically both 2/77 ~ 1 and y ~ 1, we estimate typical fermion and 
boson momenta |p| ~ \q\ ~ \J m fT <C rrif. We used that |pj <C mj, |g| <C mj 
and we have put e^ +( j ~ mj, everywhere except the exponent. In the exponent 
we used the expansion e^ + ^ ~ e p +P?Ap + <? 2 /(2e P )- Also we dropped the 
term cx 7 entering the S 10 operator, since it does not contribute to the fermion 
density and to other relevant quantities, as it is seen after the corresponding 
angular integration. 

With the help of eqs (16) and (A. 16), and also (18), (43), we obtain two 
contributions to the 3-momentum fermion distribution: 



- (Cp+ 2 ^ /70) nBoi,( ± )(p) + ^(^Hoi,^, (44) 



rf(T» + 1) / S ^Me-^^-^ S rf( To + l),„(x), (45) 
j yr/ z 



cf. eqs (22), (24), (25). Typical values of z = {t p —po)/m b in I 0s are determined 
by an estimate y 2 /2 — yr] ~ 1. The characteristic averaged value of |pj is 
~ ^j2mfT. Thereby, r] ~ 1 and thus y ~ 1. For \p\ 3> ^2rrifT typically 

y ~ 77 >■ 1. In both cases |p| ~ \q\ ~ ^ (po — e P ) 2 — mf . Cutting off the 
integral at given value of z we evaluate 



Ios(x)~-— 1 + H^/W^l+x). (46) 

X 



Here x ~ (1 + m/T/m 2 ) 1 / 2 for |p| < J2m^T of our interest, ios(^ 



23/2 

3 



X — 



1) 3/2 — >■ for x — > 1 and I 0s (x) ~ ln(x) for x > 1. 



The ratio of the fermion/antifermion density to the corresponding density 
calculated with the Boltzmann distribution, nBoi,(±)(p) = exp[— (e p ± /i/)/T], 
is as follows 



= 1 + fV* — . ( 47 ) 

PBoi,(+) Pbo1,(-) 47r z Vm fe / 

where we introduced the density of the ideal relativistic non-degenerate (Boltz- 
mann) gas 



/cftp 
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fm f T\ 3 / 2 



= pgSrexp(±^/T), (48) 

Nf = 2sf + 1 is the degeneracy factor, Nf = 2 for Sf = 1/2 fermions. 

Now we may make an attempt to solve eq. (28) in general case. We will use 
that po is near m/, and go, \o\, \p\ <^.m,f. Introducing convenient variables we 
present 
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& 2 (7o + l) 



47T 



exp(-|w|/T)0(cj), u = -p + e p , 



(49) 



where 0(u;) satisfies the integral equation: 



(50) 



Tn.iu.{uj—m b ,m b y/l+mfT/ml 



h s {u)) = <d{u - 2m b ) 



I 

mb 



de 



(f){e)J(u-e) 2 -ml 



e 2 



with e = u — u q . In (50) we separated the term leading to eq. (43) and the 
residual term. We used that e~ y / 2 (y?7) _1 sinh (yrf) ~ 1 for typical energies 
and momenta of our interest (related to y{oj) < 1) and we also cut off the 
integration in e using that y(e) Z 1 for typical \p\ Z y2m/T. For |pj ^> y / 2m^T 

one should replace ^/l + rrifT/ml in expression for I ls by ^/l + p 2 /ml. 

We may try to solve eq. (50) iteratively. First term in (50) yields eq. (43). 
Replacing this term into I ls we obtain next term of 0, etc. 

For y{uj) ^ 1 we find uj — m b <C m b and 



(j) ~ \joj 2 — ml, for T<^ml/m f , 
and for uj ^> m b , 



(51) 



4> ~ cu 



1 + 



9\ , I ■ 
-In mm 



4tt 2 



a; 



iTmf 



mi 



for T>m 2 /m / . 



(52) 



Provided g 2 /(4ii 2 ) > 1, corrections due to multiparticle rescatterings of the 
fermion are substantial already for sufficiently low temperature min(m b , T hU ) ^> 
T > ml/mf. Then, one should go beyond the iterative procedure in order to 
get an appropriate quantitative result. Thus quasiparticle approximation may 
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fail already for rather small temperatures, if the case of a warm hadron liquid 
is realized. 



For the Na interaction we estimate m N ~ g a fw, fw — 93 MeV, g a ~ 10, m a ~ 
(500 4- 600) MeV. Therefore the limiting case T <C mfjrrif is always realized 
for relevant low temperatures T <C T hU ~ m n <C m <r- (We will further argue 
that the value Tbi.f is close to the value of the pion mass ~ 140 MeV.) Then 
we may use the quasiparticle estimation of the nucleon width. For the NNa 
system for zero total baryon number ("sym"), for typical thermal momenta 
and at T ~ m 7r /2, we estimate 



n 



s P ym (p) = < ym (p) = <">1 = nJT(p) - (1-07 - 1.1)<7, 



rip yrn (p), n^ m (p), n S p m {p) and n^ m (p) are proton, neutron, antiproton and 
antineutron 3-momentum distributions. The same estimate is valid also for 
particle densities: 



p sy m=p sy m = p sym = p sym „ (L0 7 - 1 . 1 )p^. (53) 

We see that particle 3-momentum distributions and densities are enhanced up 
to ~ 1.1 times for T ~ m w /2 compared to the standard Boltzmann particle 
distribution and the density. 

Concluding, in this sub-section we have demonstrated that already at low tem- 
peratures the heavy fermion 3-momentum distribution is enhanced compared to 
the ordinary Boltzmann distribution. This is the consequence of rescatterings 
of the fermion on virtual bosons (cf. diagram (6) and its reduction to (31)). It 
results in a partial blurring of the gap between fermion- antifermion continua, 
if the case of "a warm hadron liquid" is realized. 



4-2 High temperature limit 

4-2.1 Analytic solution for fermion Green functions in STL approximation 



Assuming that typical values of \p\ ~ \J m fT are rather small (^C rrif) and 
to avoid more cumbersome expressions we further drop the term G^pj in eq. 
(35). Using (35) and (A. 2) we may as follows rewrite the Dyson equation (32) 
for the fermion sub-system derived in the STL approximation: 



G* = G°/ R + J S G°/ R \(Gf) 2 p 2 + {G R fm) + 2Gf G R mj] , (54) 
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= G)' K + J S G)' K I (Gf ) V + (G^fmJ + 2G?Gft> : 



(55) 



We introduced the quantity J, related to the operator J s from eq. (29) as 
J s . Gf(p) = J 8 Gf{p), that yields 



Js = g 2 s 



d 3 q f dqo 



(2tt) 3 7 2tt 



A b ,(+)(q)n b , {+ )(q ) + A 6i( _ ) (g)n 6) (_)(go) 



(56) 



in complete analogy with eq. (37). As it follows from (A. 5), the scalar boson 
spectral function is 



A h 



r 



(±) 



b,(±) 



(57) 



^b,{±) = ^b,(±), ReEb,(±) = ReEb,(±). For scalar neutral bosons A Sj ( + )(g) 
A s,(-)(q) and n a ,(+)(g ) = ™ s ,(-)(?o)- 

Resolving (54), (55) we find 



G 1 = G 2 (l + 2G 2 J s )- i , 



(58) 



AJ 2 S G 4 2 + 8 J s Gl + [5 

„2 ™2 



2 ~"~ ™2 



rri) 



p — mj 



Go 



m 



J s m 2 



G\ 



0. 



(59) 



Analytic solution of the fourth power eq. (59) for G<i looks cumbersome. To 
simplify the consideration we drop G\ and G\ terms in (59) (accuracy of this 
approximation is discussed below) and find the corresponding solution: 



G 



R 



fr — mj 



4J S ± 



4J S 



f 



2m). J x ( :> 



4J S 



.(60) 



This equation has the pole-like solution only for (p 2 — m 2 ) 2 ^> 16m 2 J s . It is 
obtained by taking the corresponding branch (taking negative sign in front of 
the square root in (60)) and expanding the square root term in the parameter 



C = 16m 2 f J s /(p 2 - mj) 2 . 



(61) 
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In the framework of the quasiparticle approximation typical energies are p — 
rrif ~ T and ( = 4J S /T 2 . Thus, the quasiparticle approximation works only 
for J s <C T 2 . In the leading order in ( we obtain G^ ~ (p 2 — m 2 ) -1 = G$. In 
next order eq. (60) yields 



d ~ G 2 ~ G + G 4m 2 f J s G 2 (62) 

in complete agreement with eq. (54). The proper value of the ImG^ is recovered 
with the help of the standard replacement m 2 — > m 2 — iO. In the limit case 
( > 1, G2 is already completely regular function. 

From (60) we obtain that for 12m 2 J s + 4 J s p 2 — (p 2 — m 2 ) 2 > 0: 



v 'l2m 2 J s + 4J s p 2 -(p 2 -m 2 y 
2mU(5-4 + ^) 



= - V -fr— „ 2 , 7TT ■ (63) 



For \2m 2 J s + 4J s p 2 — (p 2 — m 2 ) 2 < one gets ImGf — > 0. In the energy- 
momentum region, where the pole solution is absent, the condition ImG^ 7^ 
determines those energies and momenta, which contribute to different fermion 
characteristics, e.g., to the fermion 3-momentum distribution. 



4-2.2 Intensity of multiple scattering 

Now we may evaluate the intensity of multiple scattering J s . Supposing that 
scalar bosons are good quasiparticles in the relevant energy-momentum region, 
from (56) we find 



s ~ 9s J (2tt)3 J eio/T _ j 



(64) 



q 2 d\q\ 



2vr 2 J [ m f(T) + MT* exp (mf(T) + (3 s q 2 f 2 /T 



- 1 



where in the second line we adopted the simple form of the boson spectrum 



u 2 s (q, T) ~ mf{T) + (3 s {T)q 2 + 0(g 4 ) (65) 

which we reproduce below, see eq. (100). In a wide region of temperatures of 
our interest (5 S proves to be rather close to unity. 

In the limiting case of a high temperature typical values of momenta are 
\q\ ~ T > ml(T) and we obtain 
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Fig. 1. fo(z)/z 2 , cf. eq. (68). Dash line demonstrates asymptotic behavior for z S> 1. 

g 2 s T 2 



J, 



12(3. 



3/2 



for T>m^(T). 



(66) 



Although the STL approximation is valid only in the high temperature limit, 
let us present also the estimate of J s in the limit T -C m£ (T) in order to show 
how much this quantity is then suppressed, 



J, 



g 2 s T*/ 2 ^t b 



3/2 



23/2^3/2^ 



exp (—m* b /T) , for T< m* b {T). 



(67) 



Numerical evaluation of the integral (64) is demonstrated in Fig. 1, 



12(3. 



3/2 



TC 2 Z 2 



T 

mi 



(6S 



The horizontal line shows the asymptotic behavior (66), r s (z) — > 1 for z — > oo. 
Below we will use a simplified eq. (66) to do analytic estimates for T 3> m b (T). 
In this case the quasiparticle approximation (valid for £ <C 1) would work only 
for g 2 <C 3/3f/ 2 . The latter inequality is not fulfilled for the a meson (g s ~ 10). 

Following (36), (68) we may evaluate the density of virtual (off-shell) bosons 
in the system, cf. eq. (31), 



p b = 2ml{T){<pT) 2 = ml{T)T 2 r s /{Q(3l' 2 ). 



(69) 



For T ~ ml(T) ~ m^, (3 S ~ 1, using numerical value of J s shown in Fig. 1 
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we obtain p& ~ O.lpo- It is worthwhile to notice that virtual bosons contribute 
to thermodynamic quantities, like phonons and diffusion modes do in the 
condensed matter physics. The calculation of their contribution to the energy, 
pressure, entropy, etc is however rather non-trivial task, see corresponding 
expressions for thermodynamic quantities in Appendix B. 

4-2.3 Non-relativistic fermion distributions and density of fermion-antifermion 
pairs. 

Assume g s Z 1, T > m* b (T) and 4J S <C m 2 . Then typical 4- momenta of fermions 
are 



p 2 ~m)±0 \yfj a mfj , (70) 

as we will show it below, see after (73), (79). With the condition (70) fulfilled, 

we have |po ~ m f\ ^ m f an d \p\ ~ y / 2m/T <C mj for typical energies and mo- 
menta. This means that conditions for the applicability of the non-relativistic 
approximation for fermions are satisfied. 

If the condition (70) is fulfilled, we have 2J s G<i <C 1. Then from (58) it follows 
that G\ ~ G2- Eqs (54) and (55) coincide, if one replaces p 2 ~ m 2 in squared 
brackets in (54), (55) (but not in the pole term G°j R ). The same equation also 
follows from (59), if one neglects there a small ^J 2 G\ + 8 J S G\ term, according 
to eq. (70). 

In the non-relativistic limit for fermions eq. (60) is simplified as 



p 2 - m 2 f ± J(p 2 - m 2 f + iO) 2 - 272 
Gf ~ Gf ~ f V 2 " -> ( 71 ) 

Is 

ll = 8J s m 2 f . 
We suppressed terms of the order of 0(J s /m 2 ). 

We find from (71) that the condition A^fj~ s mf <^ \p 2 — m 2 \ that would permit 
us to use the quasiparticle approximation is not satisfied for typical values 
of fermion energies and momenta given by (70). Therefore, in the energy- 
momentum range of our interest we definitely deal with off-mass shell fermions 
described by regular Green functions. 

The accuracy of the approximation G\ ~ G 2 made to get (71), cf. eq. (82) 
below, is rather appropriate for temperatures T < T bLf . The parameter of the 
expansion is 2J s Gf (T < T hU ) < \fT s jmf < -tj=, see eq. (85) below. 
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The solution of eq. (71) that yields ImGf > should be omitted as unphysical 
one. Thus, for the energy-momentum region of our interest \p 2 — m 2 \ < \j2r{ 2 
there remains only the lower sign (— ) solution: 



2 2 

p — mi 



ReGf ~ ReGi 



7? 



ImGf ~ ImG 



R 



2 7 2 - (p2 - m j)2 



0(2 7s 2 -(p 2 -mJ) : 



(72) 



As follows from (72), the fermion spectral function satisfies the full sum rule 
(A. 12), if both regions p 2 > mj and p 2 < m 2 are taken into account. Thus, 
although we did approximations, their consistency is preserved. The Green 
function (72) is the regular function, opposite to the pole solution Gj R . 

Replacing (72) into (A. 16) we find the 3- momentum fermion distribution 



2m/\/Js/ep 



%,(±)(£) 



m/7o 



-1mf\/~Jl/e v 



7Y 



l 2 s 



exp[(f + e p =F///)/T] + l' 



(73) 



where we introduced the variable £ = Pq—^ p and used that |£| < 2mf\fT s /t p <C 
e p ~ rrif. We also have dropped the term oc pj, since it does not contribute to 
the particle density P/,(±) due to the angular integration and since Tr7 7 = 0. 
Doing further the replacement £ = —2mj\fT s je p + Ty and using that e p > 
mj ^> T, we obtain 



«/,(±)(p) 



{e p m f l +70) 



fe p T^ 2 



2nmfj! /4 



e p T 



x cxp 



e p - 2m f ^/T s /e p T 
T 



(74) 



h{x) = J dye y \jy-y 2 x 1 , x = Am f ^JT s /(e p T), 



(75) 



J 2 (x)~V/*(l-£), for x « 1, 

o z 
7 2 (x) ~ ^ , for i» 1. 



(76) 
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As we have mentioned, the condition x ~ 4y / X/T <C 1 is not fulfilled within 
the high temperature limit, which we are interested here, for g s 1, cf. 
estimate of J s , eq. (66). Thereby, we show the result in this limit only for the 
completeness of the consideration. 

Replacing (74) into (A. 19) we find the fermion-antifermion density for one 
fermion species: 



Pf,f = PM+) + PM-) ~ i e MPf/ T ) + e M~Pf/ T )} 



2T 3/2 



X 



d 3 p 



(77) 



Subsequent integration yields 



PM±) 



PBol,(±) 7rJ ; 



3/4- 



T 



cxp 



2y/T s 
T 



(78) 



m 3 / 2 T 3 /4, 



exp 



-rrif ± fif + 2\fT s 
T 



(79) 



Integrating we have set in e p ~ m/ + p 2 / (2m/) in the exponential factor and 
e p ~ m/ in other values. At this instant we are able to support our above used 
estimate of typical fermion momenta given by (70). 

In the limit ^T s /T -> (i.e. for g s -> 0, cf. eq. (66)) the ratio (78) of the 
particle/antiparticle density to the density of the Boltzmann gas (given by eq. 
(48)) tends to the unity. With the growth of the parameter ^J~T S /T the ratio 
(78) monotonously increases. Thus the density of fermion-antifermion pairs is 
exponentially increased compared to the standard Boltzmann value. 



The result (79) can be interpreted with the help of two new relevant quantities 



-},( + )( T )=-/( 1 -^ 



771/ 



2 J J s < TO/, 



(80) 



and 



TO/ TO/ 



, 2JJ S <C to/. 



(81) 



These quantities have the meaning of effective fermion and antifermion masses. 
However, contrary to the usually introduced effective masses, quantities (80), 
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(81) enter only the exponent in the expression (77). We see that m*j^(T) and 
^/(_)(7 1 ) decrease with increase of the intensity of the multiple scattering J s . 
The latter value rises with the temperature, cf. eq. (66). In case of the non- 
zero baryon density the antifermion effective mass proves to be slightly higher 
than the fermion effective mass. Thus, for /if ^ the fermion mass-shell is 
blurred a bit more intensively than the antifermion mass-shell. 

Were eq. (78) correct also for sufficiently large values of J s , we could estimate 
a value 



Js = J^l(K£ rel ) ~ (m f - ^) 2 /4, (82) 

at which the effective fermion mass (80) would vanish. °' rcl would be then 
the typical temperature demonstrating a complete blurring of the gap between 
fermion and antifermion continua. Within the non-relativistic approach from 
(66) and (82) we would get 



T^T 1 * W/ - HK 1/2 (T^). (83) 

Here the correction factor r s takes into account deviation of the numerical 
value of the integral (64) from its asymptotic value (66). 

However for T > f' rcl the non-relativistic approximation for fermions is 
definitely incorrect, since the exponential factor in (79) arose from fermion 
occupations, which in any case should be less than unity, whereas for T = 
^bi'.f reI this factor has already reached the unity. In reality, as we show below, 
the non-relativistic approximation fails at still smaller temperatures. Thereby, 
we supplied here corresponding artificial values by an additional index "n.rel" . 

The absolute maximum value of the density, which could be achieved in the 
region of the validity of the non-relativistic approximation, can be estimated 
with the help of the replacement of the exponential factor in (79) by unity. 
Then we evaluate 



3/2 T 3 3 3 W/8 m 3 / 2 T 3 / 2 

max ^ J is J 



(84) 



where we have used eq. (66). For T ~ m^, g s ~ 10, f3 s ~ 1, m/ ~ 6.7m,,- we 
estimate 2p™|- x ~ p . 
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4-2.4 Relativistic fermion distributions and density of fermion-antifermion 
pairs. Blurred hadron continuum 

An exponential smallness of fermion 3-momentum distributions disappears, if 
typical energies satisfy the condition p ± /if < T <C m,f. Thereby, and since 
we consider \/j,f \ < T, to keep the /iy dependence in equations becomes even 
less important in this energy regime. Thus we further consider the case of the 
hadron vacuum. ImG^ is non-zero at po near only for J s > mj/12, as it 
follows from eq. (63). This means that such small energies are present with 
a high probability only for J s > rnj/12. For p <C rrif the non-relativistic ap- 
proximation, which we used above dealing with the fermion energies po near 
the mass shell, becomes invalid. Thus, for J s > mj/12 we should adopt the 
fully relativistic approach in order to incorporate the region of small fermion 
energies. If po,Qo <C y/Js (qo is the boson energy variable in the diagrams 
(6), (31)), one may drop the energy dependence of the fermion Green func- 
tion at all in the calculation of the fermion 3-momentum distribution, where 
typical energies po are ~ T. The typical fermion 3-momentum is there \p\ ~ 
m N ^ \q\ ~ (T, J2mlT), see eq. (89) below. Thereby the STL approximation 
continues to hold in the given regime. 

The condition of appearance of a non-trivial contribution to the imaginary 
part of G 2 for small values p < T, 



jr(T° u )^m}/12 (85) 

determines the characteristic temperature of the complete blurring of the gap 
between fermion-antifermion continua. Within the fully relativistic approach 
for fermions, from (66), (68), (85) we evaluate the typical temperature of the 
blurring of the fermion vacuum, 



T b s 1 ,^ / 9 s 3/ V 1 ^ s - 1/2 (T b \ f ), for T° u >ml(T* u ), (86) 

for the heavy fermion - scalar boson system under consideration. If one as- 
sumes r s (T{ u ) ~ r s (Tbi f' rel ), the quantity (86) is l/\/3 times smaller than the 
artificial value (83) estimated above beyond the region of the validity of the 
non-relativistic approximation for fermions. We would like to draw attention 
to the fact that the limit T£ u <^ m^(T bl f ) is never realized, since J s is expo- 
nentially suppressed in this case, cf. (67). The opposite limit T^i.f 3> m^(T^ f ) 
can be realized, only if the bare boson mass is rather small, namely T bl f > rrib, 
cf. next two subsections. Otherwise we have 

TSu - KC*bi.f)> (87) 
and the problem of the determination of T£ l f is then a coupled-channel problem 
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of a simultaneous evaluation of quantities T£ u and ml(T£ u ). 

Now we may find the fermion 3- momentum distribution for J s > J bLf ~ m 2 /12 
(i.e., for T > Tg l f ). Let us first consider only the contribution of the energy 
region p <C \fJ~ s - F° r typical values p ~ T, for the case T <C v^s; we ma y P u t 
Po = in the expression for ImG^- Then we obtain the additional contribution 
of this energy region to the 3-momentum fermion distribution 



~ , /" dp ImGf(p = 0) (p + m fy ) 
6n u±) (p> * - j ^r-^- 



and to the fermion- ant if ermion density (one species of fermion) 



(88) 



6ptM<PoST) = - J pd® 2mH(5 + 4 + 4 4) • ( 89 ) 

/ * v mi mi ' 



S "7 



VL, = ~m) - 2J S + ^ + 16J s m 2 , 

where we used eq. (63). We see that for J s > J bLf the fermion sub-system 
represents a rather dense packing of fermion- antif ermion pairs (J s > J bLf 
corresponds to the filling of a Fermi sea, p max = for J s = J bLf ). This is 
rather similar to the standard Fermi distribution at zero temperature but in 
our case fermion width effect is significant and the Fermi momentum |p max | 
has a different value. Moreover, effective fermion and antifermion Fermi seas 
exist simultaneously in our case. 

In order to come to the quadratic equation from the fourth order one (see 
eq. (59)) we assumed that 2J|G 2 | <C 1. For p ~ and, e.g., for J s ~ 2jj bLf ) 
we have |p max | — mf/y/3 and, as it follows from eq. (60), 2J|ReG 2 | — 1/3. 
For J s » m 2 f, (ReG^I ~ 1/(2 J s ) and 2J|ReG 2 | — 1- Thus, we may continue 
to use eqs (60), (63) also in relativistic energy region (even for J s ^> m 2 for 
qualitative estimates). 

For T in the vicinity of T b s lf (J s - J bLf < J bLf ) using that |p max | = for 
J s = J s bl f and expanding in (89) all quantities in small difference J s — J s bl f we 
obtain p 2 nax = 36(J S - Jf s )/7 and 



81y/ 3T 2 {J s - J bLf ) 2 

XYlyfim) 



S Pf , {±) (po $ T) = V "° 3* ®( J * ~ J s U )- (9°) 
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Notice that the total value P/,(±) does not tend to zero for J s = J^ u ~ m^/12. 
Indeed, one should still add to (90) the contribution of the energy region 
near the mass-shell, which we have estimated above within the non-relativistic 
approximation, see (79). 

Suppose g s ~ 10, as for the aNN interaction, and rrif = ~ g s f n . Assuming 
m%(T) being significantly less than T and using (66) thereby, we would get 
T£{° { ~ 93 MeV for f3 s ~ 1, see (86). From (79) for two fermion (nucleon) 
species and for T = T^f { ~ 93 MeV, we then estimate p N ^ = 2pjj(p ~ 
rrif) ~ 4pj( + )(p ~ tTif) — O.Olpo- This is a tiny quantity. It means that in 
reality we still have ml(T) ~ at such a temperature. Thus the estimate 
T^f r ~ 93 MeV is not relevant, e.g., for the Na system, since the bare mass 
of a is (500 -J- 600) MeV, i.e. much higher than 93 MeV. We introduced a 
superscript "0" to indicate this artificial feature. 

The a meson has a large bare mass. Therefore we actually deal here with a 
coupled-channel problem, see eq. (87). In presence of nucleon-antinucleon pairs 
the effective mass of the a meson decreases that permits an extra production 
of pairs. The value of the nucleon pair density p N ^, at which the nucleon 
continuum is blurred, proves to be much smaller than the density that is 
necessary to reach the deconfinement transition at such a low temperature. 
Thus we should incorporate the factor r s < 1, cf. eq. (68). We will correct above 
estimate after evaluation of the value m* (T), see below. Also, as we shall see 
below, with inclusion of correlation effects the effective coupling constant g* 
becomes smaller than g s that results in an increase of the value T£ u . 

For artificially large values J S (T) ^> m 2 using (89), (66) we estimate ~ 
3m 2 and 



rrif 

m± > " TiwF^ -T5^- (91) 

The contribution of the energy region near the mass shell evaluated above 
within the non-relativistic approximation should be omitted in this limit. 
Thereby we replaced 5p/,(±) to P/,(±)- Also one should bear in mind that in 
derivation (89) we have put p = in the estimate of ImG, although for 
J S (T) ^> m'j the whole energy region is populated. Moreover, we suppressed 
G 3 term in (35) that might be incorrect for p 2 > m 2 . Thus (91) can be con- 
sidered only as a very rough estimate. 

Even at such high temperatures (for J S (T) ^> m 2 ) we find no end point for 
the hadron world conjectured in [5] (the value dpf^±)/dT has no singularity at 
finite T in our case). Furthermore, we see that P/,(±) may even decrease with 
the temperature increase in this energy region. Maximum available density 
of fermions can be very roughly estimated equating (90) plus (79) (in the 
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latter equation we take into account that the exponential factor should not 
exceed unity, see eq. (84)), and on the other hand (91), from where we find 
J s (pf ax ) ~ 0Am 2 f , T ~ 204 MeV, and 2pff ~ 10p for g s ~ 10. Although 
estimates of many works show that the quark deconfinement transition may 
occur at much smaller temperature for given high density, all of them are done 
within simplified assumptions. E.g., one often compares pressures of the quark- 
gluon and hadron gases to conclude about the possibility of the deconfinement 
transition. We found that the hadron phase represents in reality a strongly 
correlated state, where the number of effective hadron degrees of freedom is 
dramatically increased with the temperature following the increase of J S (T). 
Thereby, one may expect a smoothening of the transition. More likely, in 
this case the system up to rather high temperatures may represent a strongly 
correlated hadron- quark- gluon state rather than the pure quark-gluon or the 
pure hadron state. 

Note that calculating S/, cf. eq. (11), we omitted terms proportional to fermion 
occupations. These terms are as small as the ratio of the contribution of quan- 
tum fluctuations to thermal fluctuations: 



JA.rifipo + qo)^ J A 



(2ff£ J ^(2tt)4 , g2 . 

j A h n b {%)l0f fA b n b (q )^' 

We also notice that the value Jj bl ' f ) evaluated within the relativistic approach 
is smaller than the quantity J^Hi that would follow from the non-relativistic 
estimation outside the region of its validity. 



4-2.5 Boson spectrum in the regime of non-relativistic fermions. Possibility 
of hot Bose condensation 

We discussed the behavior of fermion Green functions and self-energies. Now 
let us evaluate the boson self-energy and find the boson spectrum. We will 
continue to exploit the high temperature limit (T > m^(T)) using fermion 
Green functions obtained in the STL approximation. Let us also assume, as 
before, that G 3 = 0, Gf ~ G^-, and that boson 4-momenta q , \q\ are rather 
small (go, \q\ «m/, see estimation below). Then, dropping in eq. (14) a small 
term, which does not depend on thermal fermion occupations, and expanding 
(14) in g , \o\ we find 



ReSf(g) ^ -16glJ ^ J |^(p 2 + m* f )lmG«(p) 
x[BeG*(p)n M+) (p ) + tteG«(-p)n /)( _)(po)] 
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2 f d 3 p J dp cTKeG£(gj 



d 2 ReG*(p) 



d 2 ReG*(-p) 



dp^dp v 



dp^dp 



Z — n M-)(Po) 



xq»q»(p 2 + m 2 f )lmGK(p) 



d 3 p f dpt 



(2^)3 7 (2tt) 



dReG?(p) , , dReGf(-p) , 
- P; n /i(+) (p ) + £H^ n /.(-)(P»») 



dp" 



dp'- 



(93) 



In our case according to (60) ReGf (p) = ReG^(-p) and eq. (93) is still 
simplified. 

Let us first follow the approximation of non-relativistic fermions, 1\fT 8 -Cm/. 
With the help of eq. (72), doing the replacement £ = p — e p , for £ <C e p , and 
then introducing the variable £ = — 2mjv / dJA P + Ty we get 



ReSf ~ 



2 9sPfJ 



T 



2a/X 



vx 

T 



i; 1 



T 



(94) 



where we also used that \p\ ~ y / 2m/T <C m/ and we have put e p ~ m/ 
everywhere except the exponential factor exp(— e p /T). The quantity pjj is 
the fermion-antifermion density for one fermion species, given by eq. (78), the 
parameter 



2 9sPfJ 
J s m f 



(95) 



is associated with the renormalization of the boson quasiparticle wave func- 
tion, cf. eq. (65), 



(96) 



In the limit cases we get 



I s (x) ~ ^ 2 {l - ^) , for 
16 8 

T . . 3A/7T 

j 3 (a;) ~ — ^— , tor i> 1. 



(97) 



Thus 
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ReSf(g) 



9sPf,f 



a s (qo ~ 2^)' for x 



whereas in the high temperature limit under consideration 



(98) 



ReSf(g) 



2 9sPfJ . (2 1 ^ r 4^ ^ 1 

/j- +a s (g - -g ), for x ~ » 1. 



(99) 



As follows from eq. (66), the limit rr ^> 1 is indeed fulfilled, if coupling is 
strong (g s >> 1). In the approximation (A. 20), /X//T <C 1, the real part of the 
boson self-energy does not depend on pif up to 0(fi 2 /T 2 ) terms. Comparing q- 
independent and g-dependent terms in eqs (98) and (99) we see that expansions 
hold up to rather large values of q^ and q 2 : q\-,q 2 <C \J~T s mf ~ m 2 /^/l2 for 

Assuming the validity of the quasiparticle approximation for bosons (Im£f - 
0), we find the spectrum of boson excitations that takes the form (65) with 



ft = (l --a s )/(l -a,). 



(100) 



The wave function renormalization parameter a s yields corrections to the u 2 
and q 2 terms. 



The value 



<(T)~(l-a.) 



-i 



m 



2 _ IfiPfj ( 1 _ Th{x) 



(101) 



has the meaning of the squared effective boson mass. For the most interesting 
case, x >> 1, we obtain 



ml 2 (T)~m 2 b (l-a s ) 



-i(l-%^|, 4 V / X»T. 
V mWJsJ 



(102) 



Thus the effective boson mass achieves zero at some critical temperature T c b, 
being determined by the condition 



jiP) (rp 



cbJ 



m 2 b 



(103) 



Assuming that the non-relativistic approximation for fermions is fulfilled up 



to T ~ T c s b we estimate 
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PfjVM = : ±£k • ( 104 ) 



^9s 



For T cb > T blf we have « s oc mj/mj < 1 and f3 s ~ 1 + \a s . To get expression 
(104) we used eq. (66). Note that at T = T cb the correction factor r s — 1. For 
mo- ~ 600 MeV, T cb ~ m w , g s ~ 10, a s = we estimate p^Ar — 0.7p - 

For T > T cb the value m*^{T) becomes negative leading to the instability 
of the spectrum. The stability is recovered due to the appearance of the s- 
wave Bose condensation of the classical scalar field. Such a condensation can 
be called a hot Bose condensation (HBC), since it arises for the temperature 
T > T c s b , rather than for T < T cb . As the consequence of the strong boson- 
fermion-antifermion interaction, the number of fermion degrees of freedom is 
dramatically increased that, on the other hand, results in the increase of the 
boson abundance. Boson degrees of freedom feel a lack of the phase space for 
energies and momenta ~ T and a part of them is forced to occupy the coherent 
condensate state, thereby. 

Let us now show that the quasiparticle approximation for bosons, which we 
have assumed to be valid exploiting (65), is, indeed, fulfilled in a wide tem- 
perature range. For that let us evaluate ImE^. Within the validity of the 
non-relativistic approximation for fermions, one has 2\fT s <mj. Then from 
(15) and (72) we find 



x (V + m) + pq) \j 2 7s 2 - (pi - efj ' 
1 



x 



_exp[(p - H f )/T\ + 1 exp[(p + q - fi f )/T] + 1 



(105) 



We see that in the critical point of the HBC (q = 0) the squared bracketed 
term vanishes and, thereby, the boson width also vanishes. Thus in the problem 
of the determination of the critical point of the HBC one, indeed, may use the 
quasiparticle approximation for bosons. 

Now let us consider finite but rather small values of q and \q\. For q <^ 
T, \q\ <ti \JnrifT we may drop the g-dependence everywhere except particle 
occupation factors. Separating fermion particle and antiparticle contributions, 
with the help of the replacement p = e„ — 2m f^ s + Ty, we obtain 



2^T 2 [exp^/T) + exp(-/i / /T)] (>/t _ x ) 

ImE" (q) ~ 



3/2 
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h{x) = 1 e-^dy (y - y 2 /x) = e -*(l + -)+l- -, (107) 



x 2 

h — — , for x < 1, 

J 4 ~l--, for x>l. (108) 
Finally, we find 



x I," 1 (^) e"^ {e^ T - l) • 



(109) 



Comparing (109) and (94) we see that in the high temperature limit under 



consideration, for x ~ ^> 1, and for q 2 <C T^/~T S one has |ImS^| <C 



T 

|Re£f |. The inequality |ImS^| <C |Re£f | also holds for q <C \fT s at ^r 5 - <C 1. 
Thus, £/ie quasiparticle approximation is, indeed, valid for bosons in a wide 
temperature range of our interest. For typical values q ~ T we have |ImSf | -C 
|ReS^| in both limits x ^> 1 and x < 1. Thereby, we justified that above we 
correctly used the quasiparticle approximation to calculate J s . On the other 
hand, the quasiparticle approximation fails for |ImSf (go = m t(T))\ > m£ 2 (T), 
that takes place in a narrow vicinity of the HBC critical point, but not in the 
critical point itself, where ImSf = 0. 



4-2.6 Hot Bose condensation in the regime of blurred fermion continuum. 

For J s > j]? u the fermion energy region po <C m/ is permitted. Due to that 
there appears an additional term in the boson self-energy. To find it let us 
consider the limit J s - J s blf < J s blf , for J s > J blf . Then still one has G x ~ G 2 
with a reasonable accuracy. With the help of eq. (60) from (93) within the 
same set of approximations, which we have used to obtain eq. (89), we find 



6W 7r 2 T b s lf V 4 m) 4 m}/ 
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where 5pjj(po S T) — 5p/,(+) + <5p/ 5 (_), see eq. (89). Eq. (110) yields the 
correction term to the effective boson mass (102) 



2g* aPfJ 36(ln2)^ P/i/ (po<T) \ 

where we dropped a numerically small contribution to the wave function renor- 
malization from the region po < T. 

If the effective boson mass ml(T) achieves zero in the regime T > T^ lf , T — 
Ty_ f <C T£ h{ , then the critical point of HBC is determined by the condition 



m* 2 (T) ~ ml(l - a,)" 1 



2fe/(PD - m 7 ) _ 3Q\n(2)g 2 s 5p fJ (p < T) = ^ (m) 



where we used eqs (99), (65) and (110). 

Applying results for NNa we should replace Pf j — > 2 P f j = Pnn an d $Pf,f ~^ 
2Spf j = 5p N fi in boson self-energy terms to take into account two fermion 
species, i.e., neutrons and protons in the given case. Numerical estimate shows 
that for m a ~ 600 MeV and g s ~ 10, the critical temperature of the HBC is 
T c s b ~ 114 MeV, j£ 6) ( T cb) - 1-5 J s bLf (r s = 1), a s (T c s b ) ~ 0.3. As we mentioned, 
if the limit T ^> m*(T) were fulfilled for T ~ T{ u , using (66) we would 
come back to the estimate T^{° { ~ mfj3^ A /g s ~ 93 MeV. However for such 
a temperature m* is still close to the bare mass and one needs to use the 
opposite limit expression (67). Thereby, the value T blf should be somewhat 
larger. Within a coupled channel estimate we find that the renormalized value 
T bu ~ m U T bu)- Finally we find that T blf ~ 105 MeV, that is only slightly less 
than T cb . Using instead of g s a smaller value of the effective coupling constant 
g* ~ 0.7^ s , as we estimate it below, we obtain higher values T£ u ~ 150 MeV 
and T cb ~ 160 MeV. For smaller value m a (e.g., for 500 MeV instead of 
600 MeV) we would get smaller values of T h \ i and T cb . With these variations 
we see that in all relevant cases T bl f and T cb remain to be somewhere in the 
vicinity of the pion mass (~ 140 MeV). 

Please notice that in order to find a possible relation between the HBC and 
the chiral symmetry restoration one would need to consider both possibilities 
in the framework of the very same model, e.g., the linear a-model, introducing 
meson self-interaction terms and the spontaneous symmetry breaking for the a 
vacuum at T = 0. As an intriguing circumstance, let us mention that replacing 
g a ~ m N / f n (for a — f n , as it follows from the a model) into (86) and (83) 
we find T£ u = f n and T^' rel = v^/^ for /5 S < 1. Namely in this range of 
temperatures one expects the chiral restoration phase transition, cf. [22]. We 
postpone a more detailed discussion of these questions to the future work. 
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For T > T cb , in the mean field approximation, the classical scalar field <\> sc is 
determined by the equation 



m 2 MF (psc + A? 0L = 0, (113) 
Wmf = ™& + ReSf (g -> 0) + 5ReSf (q -> 0), 

which has the solution 0;? c ~ — (m^p/Xf )0(— m^p) for T > T c s b . Here the 
value > is the effective boson self- interaction coupling constant, 



Urn = -Af 0i/4. (114) 

It arises since for T > T cb one needs to add the condensate dependent terms to 
the one term ^-diagram. With inclusion of this interaction excitations become 
stable in presence of the condensate. In reality one also has an extra term in 
the Lagrangian, L int = — A vac 0^ c /4, related to the vacuum boson-boson self- 
interaction. We for simplicity suppressed the latter. 

The condensate field is the static field corresponding to the absence of real 
scalar particles with zero momentum. The contribution of the condensate to 
the free energy density becomes 



5F(V, T) = -!?Mp- e{T _ T s h) _ (T _ T4) 2@ (T _ (115) 

demonstrating typical second order phase transition behavior, however for 
T > T cb rather than for T < T cb , as it would take place for the ordinary phase 
transition. 



4-2.7 Boson abundance in matter and at infinity 

The boson population in the medium is greatly enhanced with the increase of 
the temperature due to the decrease of the effective boson mass. Accordingly, 
the distribution function (A. 16) has a sharp peak at small momenta |g| 2 < 
m* b 2 (T)/p s : 

dN med (1 — a 

^ = ^^ = ;w_i ' for T<T - ( 116 ) 

with u s (q,T) given by (65). Distributions of particles at infinity might be sig- 
nificantly different from distributions inside the matter. This depends on the 
scenario of the breakup stage. If breakup were sudden, then particle distribu- 
tions at infinity would be given by [13,14] 
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dN°° Jml + q 2 dN med 

rf 3 g/(27r) 3 " u s (q,T)(l - a s ) A/(2tt) 3 ' cb ' 1 ' 

For sudden change of the system the particle momentum is conserved, whereas 
the particle energy might change. However the total energy is certainly con- 
served. Thus eq. (117) should be still supplemented by the requirement of 
the conservation of the total energy. The energy mismatch that arises at the 
breakup stage is compensated by the change of the energy of the particle 
collective flow. 

For T > T cb , in the vicinity of the critical point the stable spectrum of exci- 
tations is determined by the equation 



(1 - a s )co 2 MT) ~ m 2 MF (T) + 3Af {T)<t> 2 sc + (1 - -a s (T))q 2 , (118) 

for a s < 1, \q\ > 0. Here we used eq. (94) and added the contribution 
of the classical condensate field. The latter is found with the help of eq. 
(114), if one does there the variable replacement <p sc — > <f) sc + </)' to recover 
fluctuation-condensate coupling terms. To avoid a more cumbersome expres- 
sion we dropped here a numerically small correction term to the wave function 
renormalization from the region p Q < T. 

For the temperature in the vicinity of T cb the boson behaves as almost mass- 
less particle, oj 2 ~ q 2 (for a s <C 1). This results in the enhancement of the 
production of soft bosons and in the corresponding enhancement of the total 
particle yield compared to that would be for the originally massive bosons. 

Observation of a 5- function- like peak in the meson distribution at zero mo- 
mentum, if occurred, could be interpreted, as the fulfillment of the condition 
T b . up > T hc , where T h up is the temperature reached at the breakup stage. A 
significant enhancement of the meson distribution at small momenta \p\ < T 
can be interpreted SIS db SI gnal of the closeness of T b . up to the value T bc (for 
Tb.up < T bc ), see also [23]. 

The HBC may appear only, if T c s b is less than the critical temperature for the 
deconfinement, Td ec , since for T > Td cc there would occur a complete break- 
down of the hadron vacuum. As we estimated T cb ~ m n , and as we argued, 
at relevant rather small density the deconfinement is probably delayed up to 
a high temperature. Please notice that to simplify the consideration we disre- 
garded in our analysis the quark-gluon contribution to hadron quantities, e.g., 
an extra decrease of meson masses due to the change of the quark condensate 
with the increase of the temperature. 

Concluding, we treated the fermion-boson problem self- consistently. Fermions 
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Fig. 2. Logarithm of the ratio of the fermion density to the corresponding Boltzmann 
quantity for g s = 10 and g s = 7, and the effective scalar boson mass as functions 
of the temperature in units of ?y f . The vertical dash-line indicates the area of the 



due to rescatterings on bosons acquire broad widths and, as the reaction on 
that, bosons decrease their masses. 

The typical behavior of the logarithm of the ratio of the fermion density to the 
corresponding Boltzmann quantity as function of the temperature is shown in 
Fig. 2 for T < T c s b . Moreover Fig. 2 demonstrates the temperature dependence 
of the effective scalar boson mass. We see a huge enhancement of \n(pf / pf^oi) 
and a drastic decrease of ml/mb in the vicinity of T = T blf . At the critical 
point m* bt /m b reaches zero demonstrating possibility of the second order phase 
transition to the HBC state. 

4-3 Contribution of baryon resonances 

Above we considered a simple example of the system consisting of heavy 
fermions of one species interacting with one kind of less massive scalar bosons. 
However at finite temperature together with nucleon/antinucleon states the 
high-lying baryon/antibaryon resonances, like A isobar, hyperons, etc, are 
also populated with some probability. Resonances interact with each other 
by boson exchanges, as well as by a residual interaction (including repul- 
sive baryon-baryon correlations). To describe the multi-component system of 
the baryon/antibaryon resonances interacting with mesons we need to know 



HBC. 
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coupling constants between different particle species. In general, Dyson equa- 
tions for Green functions of different particles are coupled and the problem 
proves to be very complicated. Different meson exchanges may significantly 
contribute. E.g., for the A isobar the coupling in the pion channel is the dom- 
inant one. However the scalar channel might be also important. According 
to [24] a mesons interact with A-isobars with the same universal coupling 
constant g s ~ 10. In case of iV*(1440) one gets [25] g a NN* / QaNN — 0.47, 

9nNN* 1 9-rrNN — ^ -^daN N* / daN N — 0.25. 



4-3.1 Low temperature limit 

As above, let us for simplicity first assume that baryon resonances couple by 
an exchange of only a scalar neutral boson (s). This simplification is sufficient 
to find particle distributions in the low temperature limit, T <C min{mb, Tbi.f}. 
Then we may still use eqs (A. 16), (A. 19), (16) for the given baryon resonance, 
however with A° operators being different in dependence on the spin of the 
baryon species. Calculating the resonance width we need to take into account 
in (12) various possible intermediate states, since the given resonance may de- 
cay to the virtual boson and to an another baryon resonance then absorbing 
back the virtual (off-mass shell) boson. Notice that we discuss only temper- 
ature effects. Just to simplify the consideration we artificially suppressed the 
widths terms surviving for T = 0. 

Let us present the density of the baryon resonance/antiresonance states of the 
fixed species B*. In the low temperature limit using (20), (47), (48) we obtain 



Pb*,{±) — N B * 



m B * 
2tt 



jk 3/2 



exp 



m B * =F Hb* 



T 



(119) 



fm B *T\ 3/2 



m B * =F Hbt ' 


9sB*B* 




I T 


to* l0S 


V m b J 



N B * is the degeneracy factor (e.g., N n — N p — 2). The summation is performed 
over all possible states including the given baryon state (B*), m B * is the mass 
of the B* baryon resonance. For non-strange baryons with the same baryon 
number, as for the nucleon, we have fj, B * = /in, if there is a permitted reaction 
channel: B* <-> N + b mrt . Here b mrt is the virtual boson (in our model example 
it is a scalar boson, whereas in reality it also could be the pion, a, u, etc). 
Second term in (119) is due to the diagram 




where external double-lines correspond to the given resonance B* and internal 
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solid line relates to the resonance B*, including the given resonance B* and 
the nucleon state. In the latter case the self-energy term is the same as in eq. 
(6). 

Thus, the baryon resonance density is substantially increased compared to 
the Boltzmann value. The term in eq. (119) relating to the decay of the given 
baryon resonance to the nucleon and the virtual boson (i = n or p, solid line 
in (120)), has no suppression factor, as e~^ rnB *~ mN ^ T . The latter factor would 
arise, if one worked in the framework of the quasiparticle picture. If N B *g 2 sB * N 
were larger than Niyg 2 NN , the density of a high-lying resonance could be even 
higher than the nucleon density, showing a laser effect. Opposite, if g S B*B* were 

1 3 

negligible for % ^ j, the laser enhancement would disappear and a high-lying 
resonance state would be less populated than the nucleon one. Nevertheless, 
in any case the resonance state proves to be more populated compared to the 
value determined by the corresponding Boltzmann expression. 

4-3.2 High temperature limit 

To proceed in the high temperature limit let us additionally assume that 
9sB*b* 7^ only for % — j. Then in the non-relativistic approximation for the 

% 3 

resonance, the density of a B* baryon resonance (and its anti-partner) is found 
with the help of eq. (79). We obtain 



" vi^r xp i — —\ h {—)- (m) 

Here, in agreement with eq. (56) the intensity of the multiple scattering is 



^ ' o 

where we used that for scalar neutral bosons A s ^ + )(q) = A s ^(q) and n S) ( + )(g ) = 
n s (_)((?o). The effective mass of the baryon resonance follows from expression 
(80): 



m* B , ( ±)(T) = m B * 1 =F — : — , for 2JJ s B * < m B *. (123) 

' v ' I m B * m B * J v 

Analogously, one rewrites expressions (90), (91) for the density of the baryon 
resonance in the relativistic energy region. Within the quasiparticle approxi- 
mation for the boson using (99), (110) we obtain the effective boson mass: 
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mf (T) ~ (1 - a,)' 1 (124) 



x 



m. 



v ,'' 2 9sB*b*Pb* (Po ~ rn B *) 36{ln2)g 2 sB » B *5p B: (p < T) 

E — 1 • A— + 



s 

bl.f 



E 7 ' -» (125) 

J s ,B*mB* 



The summation is over all baryons and antibaryons, p#* is the corresponding 
baryon or antibaryon density. As in eq. (Ill), we dropped a numerically small 
correction term to the wave function renormalization from the region po < T. 

In the limiting cases of high and low temperatures we recover eqs (66), (67), 
now with the coupling constant g a B*B* standing instead of g s and 0b*, instead 
of ft, 



Js,b* = , for T»m*(T), (126) 



and 



2 T 3/ 2 . , 

^,b> = 2 l/ v/2 ^2 b - eM-m* b /T), for T«m*(T). (127) 

We see that for T 3> rn* h (T) the intensity of the multiple scattering, J Sj _b*, of a 
high-lying baryon resonance would exceed that for the nucleon, J Sj at, if g 2 sB *B* 
were larger than g 2 NN . As follows from eq. (123), for g 2 B *B* / m B* > 9sNNl m2 Ni 
the baryon continuum for the given high-lying baryon resonance would be 
blurred at a smaller temperature than for the nucleon. Please notice that 
such a relation between coupling constants is not fulfilled for the realistic 
hyperon-a-nucleon interaction, cf. [27]. Nevertheless, using estimates [27] we 
conclude that resonances contribute essentially to the total baryon-antibaryon 
density for T ~ Tu.f. Moreover, we artificially suppressed all couplings except 
9sB*b* for i — j that is certainly not the case in the reality. Also in reality 
9sB* b* 7^ 9sB* b* for % 7^ j since N B * ^ N B * , that may stimulate in some cases 
a laser effect. A high-lying state might be more populated than a low-lying 
state (Note that, e.g., in case of the A isobar, i.e. (|, |) resonance, one has 
Na/N n = 4, that works in favor of the laser effect). 

Concluding, indeed, we deal with the hadron resonance porridge at a suffi- 
ciently large temperature. 
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4-4 Evaluation of correlation effects 



Above we discussed properties of the system described within the simplest $- 
derivable approximation with only one diagram (2). Exact fermion and boson 
self-energies are determined by diagrams (4) and (5) with one free and one 
exact vertices. In the low temperature limit vertex corrections are negligible. 
Thereby, we further consider the high temperature limit T > m^(T). The equa- 
tion for the vertex can be greatly simplified within the ladder re-summation: 




Jjr^ = \ — + 

/ / ^ (128) 

that reads as 

Mad(<?i,P,P + gi) = V (qi) - J G f (p- q)G f (p + q 1 -q) 

xiV\ad(q,p + qi-q,p + qi)G b (q)v lad (-q,p,p - q) 

d 4 q 

xViad(gi,p- q,P + qi ~ ( 129 ) 

the corresponding matrix indices are implied. We need — , — , — and +, +, + 
vertex functions of the same signs, since in the formalism that uses full Green 
functions, cf. [21], any extra full Green function G h or G + ~ corresponds to 
the scattering process involving extra particle in the initial and the final state. 
These processes are suppressed, if the number of fermion- ant if ermion pairs is 
not too large. In the STL approximation, which we now exploit, we find 



~ 1 + Js(V Ud -)(P 2 + + 2M) [ReGxG*)] 2 ' ^ 

where we also used the spin structure of fermion Green functions (35), within 
the ansatz G 3 = 0, assuming Gi ~ G 2 - The quantity J s (V[~d ) ^ s Js with 
the bare vertices g s replaced to the full vertices, see corresponding fat dots in 
(128). Thus, one may restrict the consideration to the first diagram (2) only 
if (p 2 + mj) J S (V^-) [ReG 1 (p)] 2 « 1. 

At T < T^' f using (60) we evaluate (p 2 + m))J s [ReGi(p)] 2 < 1/2. For artifi- 
cially large temperatures T>Ty f we would have (p 2 + m 2 )J s [ReGi(p)] 2 ~ 
"/? (2./J < 1. 

In the full series of vertex diagrams, beyond the ladder approximation, there 
are graphs with crossed boson lines. In the STL approximation each diagram 
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that includes the same number of full scalar boson Green functions yields the 
very same contribution independently on where the boson lines are placed 
inside the diagram. Counting the number of boson lines with full vertices in 
first diagrams of S j we find 



1 + J S (V )(p 2 + m} + 2pm f ) [Red (p)] 2 

| g a J!(V— )(P 2 + m} + 2j>m f f [ReGi(p)] 4 
1 + QJ s (V—)(p 2 + mj + 2pm f ) [ReGi(p)] 2 ' 

Thus, the ladder approximation yields an appropriate estimate of the full 
vertex up to rather high temperatures. 

For rough estimates we, as before, may consider only one diagram of $ but 
with an effective coupling constant g* instead of the bare vertex g s . For low 
temperatures we have g* — g s . The vertex suppression factor increases with 
the temperature reaching the value g* ~ 0.7 g s for T = T£ u and again g* — > g s 
for T 3> Ty f . Note that the vertex suppression factor essentially depends 
on the structure of the fermion - boson interaction. For the Nir interaction 
the corresponding vertex would be less suppressed, cf. [5] and a discussion in 
subsection 6.4. 



Note that we discussed just a model example. We suppressed a possible boson- 
boson self-interaction. Inclusion of the latter complicates the consideration 
yielding a repulsion [28,26,14,29]. Also simplifying we considered only the 
baryon interaction with the scalar meson disregarding the baryon interaction 
with other meson species. A residual baryon-baryon interaction has been also 
dropped out as well as an interaction between different meson species. Fur- 
ther we will proceed step by step studying relevant models, thus permitting 
different types of interactions. 



5 System of heavy fermions and less massive vector bosons. 



Now we will consider another example, the fermion - vector boson system 
with the coupling given by 



J int 



-g^W- ( 132 ) 



The bare vertex is Vq = gv^f^. Then we will apply results to the NNu and NNp 
systems. Again we first solve a model problem disregarding other relevant 
couplings of vector mesons with other mesons, tensor coupling with nucleons, 
etc. 
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5. 1 Spin structure of vector boson propagator 



In the medium the Green function of the vector boson changes as follows 



L tl1 j T ltlJ hnh^ 



{Gb V 'CP-mi- Si + a 2 - ml - El ml ' 



i" = (g-u)g p -uV 1 h fX = q»/ s Jq~2, q 2 > 0, 

■u is, as before, the 4- velocity of the frame. In the rest frame one has u = (1, 0). 
The retarded self-energy of the vector boson is subdivided to the longitudinal 
(Ej,) and the transversal (E£) parts 

££" = ££L"" + E* 6 T^. (134) 

One usually assumes h li E^ 1/ = 0. For Ej,, E£ — > 0, eq. (133) coincides with 
(A.4). 

From (133) we find 



= -£(E?)8, = ^(Sf )J, <,J = 1,2,3, (135) 



£(E b R )*(g = 0, |g] - 0) = 2(Ef ) t (g = 0, \q\ - 0), i = 1, 2, 3, 

i 

(E?)8(«> = 0, |9l - 0) = (Sf )'( 9o = 0, |g] - 0). (136) 
5.2 Low temperature limit 



The low temperature limit, T <C min{m b ,T W j}, is considered quite similar to 
that for scalar bosons. We replace in (20) 



Vo(q)A f (p + q)V (-q)A b (q) 
9v\ 2 P + 2 ®- 4m f + - — 2 + — j + -V • (137) 



m b m b m b m b 



Assuming that p ~ m,f, q 2 = q 2 — q 2 ~ m\ and dropping linear terms in 7 
and the term oc pq, which do not contribute, we get 
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S° * 3^(70 - 1) + 9 2 V ^P^ + 3A7o. 
Then from (20) we find 



'13* 



rf(po - e p ) ^ gexp - Po )/T] ^^e- 2 / 2 ^ 



e„ r - mi 



yr) 



x (3(70 - 1) + %(e P - Po) 2 - m 2 b ] + 3 7o (e P - Po)™/ } , (139) 
I m b ) 

for e p > po, y and 77 are the same as in (43). 

Using (139) and (24) the fermion distribution is presented as follows 



(e p + 772/70) /-k 
n/,(±)(p) - TT- 1 ™Bol,(±)(pJ 



8tt 2 



™Boi,(±)(p), 



(140) 



J 1/7/ Z 2 



(141) 



4(^)-/ Sinh J y?7) e-^ 2 (. 2 -l) 1/2 ^. 



(142) 



Here variables y and 7/ are determined as in (25). Cutting off integrals at z 
corresponding to y ~ yrj ~ 1 we estimate 



(x 2 - I) 3 / 2 3 
— h 

x 2 



2 /2 31nQr+(x 2 -l)V 2 ) 



j£,(a;) ~ (x 2 - l) 1 / 2 + arcsini - 

x 2 



(144) 



with x ~ (1 + m/T/m 2 ) 1 / 2 the same, as in (46), and JoiX 2 - ~~ * 00 ) ~~ * a;2 /2, 
- 1) - 2^(x - I) 5 / 2 , - 00) - x, - 1) - - I) 3 / 2 . 

Finally, the fermion particle and antiparticle densities are 
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2 



PBol,(±) 



~ 1 + 



4tt 2 



7° 

j 0d 



\to&/ 2m f rrib 



(145) 



Correlation effects are negligible in the low temperature limit. For T <C to|/to./ 
the term oc is the dominating term. Compared to the scalar boson case, 
see eq. (47), here boson distributions are 3^mb/(2^m/) times suppressed (for 

9v ~ 0*, ™& < and ~ 1 + . For T > to^/to, the term 



PBol,(±) 

oc /£, is the dominating term and P/,(±) ~ 1 + 9 j"^ f 1 . Thus the enhancement 

^ to PBoi,(±) 87r 2 m^ 

is here higher than in the scalar boson case (again for g v ~ g s and rrif m h ). 

One can easily show that, as for scalar bosons, in the vector boson case the 
quasiparticle approximation fails for the description of the warm hadron liquid 
of a small fermion chemical potential for the relevant value g v ~ 10 ; to/ >> to&. 

For the NNu system of zero total baryon number we have to w ~ 782 MeV. 
The g^ coupling is less known. Its evaluation used in the relativistic mean field 
models [30] yields g u ~ 8 4- 10. The mass of uj is rather high. Thereby, the 
limit To;, T ml/rrif is not realized. For T ~ to^/2 (such a temperature 
is in the range of a slightly heated hadron liquid for uj) and for g u ~ 10 we 
estimate p^ = pT = pT = fiT * I-VbT- 



5.2.1 Vector-is o spin vector boson - fermion system 

For the vector-isospin-vector boson - fermion coupling (Np sub-system) the 
interaction term of the Lagrangian density is given by 



— — * 

Lint = -Qi.v.^l^r^i). (146) 

We suppress possibility of a tensor coupling. The latter has been discussed 
in [31], where the nucleon-antinucleon loop diagram has been studied for the 
case of cold nuclear matter. Quite similar to the vector boson case we obtain 



n f ,(±){p) ^ - n BoU±) {p) 



N- a 2 
8tt 2 



2e p 

T° < 6p \ _l_ T 1 ( 6p \ 

Iov[ m b ) + 2m/ 0v[ m b ) 



(147) 



WBol,(±)(p)- 



In comparison with eq. (140), here appeared extra isospin vector degeneracy 



factor Ni_ v = 3. 



The fermion particle and antiparticle densities are 
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PU±) _ n , N i-v9: 



PBol,(±) 



~ 1 + 



2 



4tt 2 



rO ( m f \ , 3m b rl /™A 

0u \m b ) + 2m f 0u{ m b [ 



(148) 



The p meson mass is 769 MeV and the coupling constant is g p ~ 3, cf. [31]. 
With these values we obtain one fourth of the u meson density, for g u = 10. 
Summing up contributions of a, u and p at g s — g v — 10, gi, v = 3, T ~ m n /2, 
m a ~ 600 MeV, we arrive at the estimation p s p ym ~ p^ m ~ p s p ym ~ p s ^ m ~ 



5.3 High temperature limit 



5. 3. 1 Analytic solution for fermion Green function in STL approximation 

For high temperatures, T > m* bt {T), where m^(T) is the effective mass for 
the transversal vector boson, we may use the same STL approximation, as 
we exploited for the scalar boson. In the latter case we first considered the 
problem within one diagram of $, i.e. without inclusion of correlations, and 
then estimated the contribution of correlation diagrams. For vector bosons 
within the STL approximation we may solve the problem in general case using 
the Ward - Takahashi identity 



QQ- 1 

-a^ = V»(q = 0,p,p), (149) 

OPp, 

where V th is the full vertex function. 

The full Dyson equation in the STL approximation takes the form 



8f = eo /S _ e« j r7 ,,e ? ^i_c« (150) 



or equivalently 



(^-^ = 1 + ^7,^. (151) 

Here J% v demonstrates the intensity of the multiple scattering, which is now 
the tensor function, 



JtT = 9t j ^ / ^ H( + m + )(9o)T^ + A l bM n b , i+) ( qo )L^ 
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+gl S 7^ Jfr K-^H^r+^HW^] » ( 152 ) 



(27r) 3 7 2tt 



and there appear transversal and longitudinal spectral functions 



A l(±) ~ 



-2ImEjf ±) 



g 2 - m 2 - ReS^ ±) 



ImE bf±) 



r, j = {t,/}, g >0. (153) 



The general structure of the Green function is determined by eq. (35), where 
we again assume G3 = 0. Using it we work out the tensor structure of (151): 



K(±)^(±) t ;+4(±)^(±) l m) ^ 



dp 2 



dd 
dp 2 



= 4, (±) n 6>(± )(^- W 

-(4 (±) - 4, ( ±))^(±)^^ 



^ 2 + -^-m f 2p u + Gi7„ 



dd . dG 2 „ ' 

-^2p4+—m f 2p u + G llv 



Here 



(154) 



dd 



dd 



dp 22pJ + lp^ m ^ + Gl ^ 



dd 



Qri , of 2 iHM)^dG l / j(j>g) \ 



(155) 



2^ + -—rm f 2p u + Gi7„ 



dp 2 dp 2 



= G 1 



(156) 



+2 



(gu ■ qp — up ■ q 2 ) 
q 2 [q 2 — (qu) 2 ] 



-^(q u •ii>-ii>- q 2 ) + ~ d ~^ m f^ u ■ ^ ~~ ^ 2 ) 



To avoid more cumbersome expressions we used a symbolic notation 



9 dGi, 2 



dp" 



2 , whereas in general case G 12 depend separately on p and p. Approxi- 
mate equality holds only for non-relativistic fermions and for T 3> m^ 2 (T)/m/. 



In the rest frame for the hadron vacuum case (/ij = 0), dropping the term oc 7 
and the linear term oc qp, which do not contribute to the particle densities, 
we rewrite eq. (156) as follows 
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Gi-2 



+ ^r' 7 ° 



(157) 



To simplify expressions we also suppressed contributions oc p 2 dG , x i which are 

small for T 3> rn* b 2 (T)/rrif. In these assumptions from (151) we derive two 
coupled equations 



dG 

p 2 G 1 - m 2 G 2 ~ 1 + (3J * + J^Jd - 2( Jl + Jr^P^, 



dG 

G 2 -G 1 = -2(^ + ^ 1 )-4. 



(158) 



We introduced 4 types of intensities of the multiple scattering 



<i 3 g f dq, 



(2tt) 3 y 2tt 



4,(+)™6,(+)(<?o) + 4,(-)^,(-)(<?o) 



, / g 2 d 3 q f dq 
J i — 9 V 



T" - a 2 
J oi — 9v 



q 2 (2tt) 3 J 2ir 
d 3 q °f dq 



( A t(+) n b,(+)(qo) + K,(-) n b,(-)(Qo) 



(2tt) 3 7 2tt 



x 



J ll — Uv 



«(+) - 4,(+)K(+)(<?o) + «(-) - 4,(-)K(-)(<?o) 
q 2 d 3 q f dq 



x 



g 2 (2tt) 3 7 2tt 
«(+) - 4,(+)K(+)(9o) + (4,(-) - 4,(-)K(-)(?o) 



(159) 



Compare these expressions with eq. (56). The consideration is simplified, if 
one assumes ^,(±) — ^6,(±)' being correct in case of a small change of the 
effective boson mass. Then and can be dropped out. 

Let us assume the validity of the quasiparticle approximation for vector bosons 
and use the spectrum 



u/(q, T) = m* 2 (T) + (3 v {T)q 2 + 0(g 4 ) 



(160) 



for transversal modes with m^ 2 (T) > 0, f3 v > 0, cf. (65). Then in the limiting 
cases T < m* bt (T) and T > m* bt (T) we find 
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Fig. 3. f\(z)/z 4 , cf. eq. (163). Dash line demonstrates asymptotic behavior for z S> 1. 



30{3 5 v /2 mlf 



for T»m*,(T), 



(161) 



3g 2 T5/ 2 



^exp(-m* bt /T), for T «C m* bt (T). 



Numerical evaluation of the intensity is demonstrated in Fig. 3, 



j. = W ^,), r .w ^ 15/1(2 



7?) 



(162) 



(163) 



The horizontal line shows the asymptotic behavior (161), r v — > 1 for 2 — > 00. 
Comparing (161), (162) with expressions (66), (67) for the scalar boson case 
we see that in the high temperature limit (T ^> m^(T)) the intensity of the 
multiple scattering J" is much higher in the vector case (for g v ~ g s and 
for m* u ~ m*). At these conditions the quantity Jg is of the same order of 
magnitude as Jg. Please notice that — > 00 for m^(T) — > 0. Thus we deal 
here with a strong critical opalescence. 

Replacing G\ from the second equation (158) to the first one we obtain one 
differential equation 



(pg -el- sr - r m )G 2 - 1 + 2{ji + r u )(2 P 2 - sr - r 01J 



dpi 
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+4(J[ + J n) 2 Plj^ = 0- (164) 

Equation (164) can be solved perturbatively only for very small values of 
Jq, J", Jq 1? J^. For typical values pi — e 2 <C mfT permitted within the 
quasiparticle approximation, the perturbative regime would hold only for 
J\,Jqi,J\i <C T 2 . On the other hand, as we have mentioned, at rather low 
temperature the STL approximation is not applicable anymore. Even in the 
low temperature limit for the regime of a warm hadron liquid the quasiparticle 
approximation (for g v ~ 10) and the perturbative approach fail, as we have 
shown it above. 

5.3.2 Non-relativistic fermion distributions and blurring of the fermion con- 
tinuum 

In non-relativistic approximation for fermions, as we shall see it below (cf. eqs 
(167), (171)), typical fermion 4-momenta of our interest are given by 

p 2 ~ m) + O (yj?m f ; J^mfT' 1 ^ . (165) 

For such momenta one may drop the higher order derivative term in (164) 
Then eq. (164) is simplified as 



^Y^^ + ^ 2 -l = 0, zcpl-el, (166) 
( 7 - c1 ) 2 ~ 8(J[ + r n )m}. 

The form of eq. (166) coincides with that have been used for the description of 
the blurring of the electron gap in semiconductors due to the electron-phonon 
interaction [6]. Heavy fermions (e.g., nucleons) play the same role, as electrons, 
whereas light vector bosons (e.g., u with an effective mass m* t (T) <C m N ) play 
the role of phonons. 

Solution of (166) satisfying appropriate condition G 2 — ► 1/z for z — > oo (z ^> 
7^ rel to be more concrete) can be presented in the integral form 



oo 

G ^) = 7^ J(z-vre-^-»d n . (167) 

v — OO 

Now we may recover the explicit expression for the full vertex V^(q = 0,p,p). 
From (158) and (166) we find 
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1 — yC 

G f ~G^ + m f )+^-^, (168) 

i.e. Gf ~ G 2 {i> + rrif) for typical energies and momenta of our interest. Using 
(149) and (166) we obtain 



In order to calculate fermion distributions and the density we need to know 
ImG. From (167) and (168) we find 



l m G 2 = -^e- z2 l^ l )\ l mGl = l mG2 1-A> . (170) 



e ' si " ' , mil 

We see that ImGi ~ ImG2 for \z\ ^ m 2 of our interest. The imaginary part of 
the Green function given by (170) satisfies the sum rule (A. 12) for \z\ <C m 2 . 

For typical energies given by eq. (165) with the help of (A. 16) we obtain 



n U±) - 



7o + l 



rcl 



exp 



(e P =F Hf) 



7o + 1 



-cxp 



£ p =F Vf) 
T 



+ 



T 

(7°' rc1 ) 2 
16e 2 p T 2 



dzexp 



z 



,(7r e1 ) 2 



+ 



2e p T / 



(171) 



Compare this result with (74). For y/Jf + Jf[ <C T (low temperatures) eq. 
(171) yields only a slightly corrected Boltzmann distribution. However (171) 
becomes incorrect in the low temperature limit due to the violation of the 
STL approximation in this limit. In opposite limit case, yj J\ + Jfi ^ T, the 
3-momentum fermion distribution is exponentially enhanced. 



For the fermion- ant if ermion density we obtain 



JV _|_ JV 

Pf,(±) = PBoi,(±)exp ( 1 2T2 11 ) . (172) 



The value 



TV _|_ JV 

m f =m f -n f (173) 

as it is extracted from (171), plays the role of the effective fermion mass. 
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The temperature T^' re is estimated from the condition m*j(T^' TC ) = 0. Then 
from (173), (163) we evaluate 



rpV,ll.rcl 

J bl.f 



60m* bt 2 (m f - n f ) 



71 9v r v{J- bl.f ) 



1/3 



(174) 



Simplifying, we suppressed the contribution in this estimate. Comparing 
(174) with (83) for scalar bosons for r v ~ r s ~ 1 we see that in the heavy 
fermion limit (rrif ^> ml) and for g s ~ g v we have f' rcl < °' rel , i.e. in this 
limit vector bosons contribute more in the blurring of the fermion continuum 
than scalar bosons. Please notice that deriving eq. (174) we assumed A^(±\ — 
A l b ^±y The latter condition is not fulfilled for m* ht <C m b . Thus (174) can be 
considered only as a rough estimate. 

For the uj meson from (174) we would obtain T£{° f ' n ' iel ~ 2.3m,,-, for g w = 
10, r v ~ 1 and for the mass m* bt = m w = 782 MeV. But then T^{° f ' n ' rel < 
m w(^bi'.f' nrel ) an d the condition for the use of the high temperature limit is 
violated (r v 1). By index zero we again indicate that we will refuse of 
this estimate. Thus one should take into account the decrease of m*, + for T ~ 

Cut 

Ty'f' re . Taking m* ht ~ we would get " ' rc ~ 120 MeV. However, as we 
have demonstrated in case of the a meson, fully correct estimate should be 
done within self-consistent account of both the fermion and the boson coupled 
channels. Finally we may roughly evaluate T^' vel ~ m vr- 

We supplied the value T^'"'™ 1 with the index "n.rel", since the estimate was 
done in the non-relativistic approximation for fermions. In pure relativistic 
case one needs to solve much more involved eq. (164). Without additional 
calculations we may only conclude that for relativistic fermions the value f 
is still less than T^'°' rel . However in our model example we did not take into 
account interaction with other mesons, e.g. scalar boson correlation terms. 
Such terms suppress vertices that may result in an increase of the actual value 
of T b \ f . 



5. 3. 3 Vector-is o spin vector boson - fermion system 

The vector-isospin vector boson - fermion system (the interaction term Li nt = 
—gtv^luT^^) is considered in the same manner, as the vector boson - fermion 
system just discussed. All expressions obtained for vector bosons continue to 
hold also for vector-isospin vector bosons, if one replaces g v — > g* ~ V3gi. v .- 

For uj and p mesons one has g w ~ 8 4- 10 and g p ~ 2.5 4- 3, from where we 
estimate that the effective p meson coupling is ~ 0.5 of that for the uj. 
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5.3.4 Comparison with example of one diagram of $ (no correlations) 

To compare our general result with the result one would obtain without in- 
clusion of correlations let us consider the problem with only one diagram of $ 
(first diagram (2)). In this case one should replace the vertex (149) by the bare 
vertex 7 M . In the STL approximation the Dyson equation (150) then takes the 
form 



Gf = G°/ R - G°/ R Jr^Gf lv Gf, (175) 

with given by eq. (152). Again, to simplify the consideration we will assume 
that A\ ~ A\ being correct in the case of a small change of the effective 
vector boson mass. This assumption allows us to omit values and J\ x . 

Now we may work out the tensor structure of (175). For that we first calculate 
the auxiliary quantity 



K - ^f)^^i> + Gzm f ) lv (Gij> + G 2 m f ) (176) 
= _ G y + 2m f p'G 1 G 2 + 3m 2 f G 2 2 - 2^^G\ - 2m f ^-G 1 G 2 , 
see eq. (155). 

Let us again use the non-relativistic approximation for fermions. Then we may 
suppress terms containing pq compared to the corresponding term oc p <?o ~ 
rrifqo. Also we drop the term oc 7, which does not contribute to the particle 
density. Then we obtain two coupled equations 



p 2 G 1 - m)G 2 



1 + J v Q {G{p 2 - 3G 2 2 mj) + 2J v 2 m)G{, 



G 2 — G\ — 1J\G\G 2 , 



(177) 



with 



fj2^ J 27 K(+) nfc '(+)^°) + At b,(-) n b,(-)(lo)} ■ ( 178 ) 

Within the quasiparticle approximation for bosons J 2 is expressed through Jg 
and J%. 



We replace G\ from the second eq. (177) to the first one and retain there 
only quadratic terms in G 2 , where we put p 2 ~ m 2 , and we use J± <C m 2 . 
Approximations are correct for typical momenta of our interest given by (165). 
Then we find 
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p 2 - ml ± i A /2( 7 ?- re1 ) 2 - (p 2 - ml + iO) 2 
G2- " V 7 lv> -■ (179) 

( 7 n.rel)2 

For \p 2 — m 2 \ < \^2^- Tcl the upper sign solution should be dropped. This 
result is nicely matched with that in the scalar case, cf. (72). The difference is 
only in the value of the intensity of the multiple scattering (J% in the vector 
boson case instead of J s in the scalar boson case). Imaginary part of the Green 
function (179) satisfies the exact sum rule. However we note that the Green 
function (179) is quite different from that given by eqs (167), (170), which 
exploit the full vertex. 

We conclude that in the vector boson case correlations qualitatively change the 
picture modifying expressions for the effective fermion mass and the value of 
the temperature T£ u . 



5.4 Hot vector boson condensation 



Now let us consider a possibility of the HBC for vector bosons. To sim- 
plify the consideration we will not discuss a modification of the vector boson 
spectrum rather we will only estimate the critical temperature for the HBC 
in the transversal channel. For that we need only to evaluate the quantity 
Re(E 6 )j( go = 0,|g1->0),cf. (136). 

Using (14), where we replace one bare vertex to the full vertex (149), we obtain 



M^riQo = 0, \g\ - 0) ~ -4^ 2 Tr / ^- J^YReGf(p) 

^ ' 



dp^ 



x Re—^lmGf(p)n f ( Po ) + {/ - /}. (180) 



With the help of eq. (169), where only second term contributes, assuming 
G\ ~ G 2 (non-relativistic nucleons) we find 



R«(fif)!(» = 0,1^0)= ^Ij^j^ 

v ' 

xlmG^po) + {/ - /}■ (181) 
Using (A. 16), (A. 19) and (136) we obtain 



Re(Sf Y(q = 0, \q\ - 0) ~ -6g 2 vPfJ /m f . (182) 
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Here pf j is, as before, the density of fermion-antifermion pairs of one fermion 
species. In reality other fermion species (e.g., as neutrons and protons) may 
also contribute to E b . 

The critical temperature and the fermion-antifermion density for the HBC of 
vector bosons in the transversal channel under consideration would be found 
from the condition 



m 



* 2 = m\ + Re(Sf ) tr (q = 0, \q\ - 0) = 0, (183) 



if it had the solution. However, since the intensity of the multiple scattering 
Ji tends to infinity for vriQ — > 0, see (161), the density of fermion-antifermion 
pairs given by eq. (172) is anomalously increased. This motivates the pos- 
sibility of the first order phase transition to the HBC state with a jump of 
m^ 2 (T c k) from a positive to a negative value in the critical point. 

Let us assume that this jump occurs at rrCQ = C m bv < C < 1- Then with 
the help of eq. (182) we find 



Note that eq. (172) was derived in the non-relativistic approximation for 
fermions. It should be replaced to the relativistic expression at temperatures, 
when m*f(T) becomes to be < T. Thus eq. (184) yields a correct estimate only, 
if m*f(T) > T. To further study the first order phase transition one needs to 
calculate the contribution to the thermodynamic potential in both phases. 
This requires a more detailed analysis. Only for a rough estimate of the value 
T c b one may use (172), (161), (184). 

If the system consists of two fermion species (e.g., neutrons and protons) 
it can be taken into account by the replacement p^y — > 'Zp'fj — P* n m 
the boson self-energy. For the NNu (two fermion species) and for g v ~ 10, 
~ 782 MeV, m N ~ 933 MeV, assuming m* bt (T^ h ) ~ we estimate 
p*^ ~ 0.7p . According to (172) and (163), the density p*^ ~ 0.7p cor- 
responds to the temperature T ~ 110 MeV. The value 0.7po at such a low 
temperature is well below the deconfinement density. Thus quark effects seem 
to be still unimportant. However note that values Ty'f' rel — 120 MeV and 
T£ h ~ 110 MeV could be underestimations, as the result of crude approxima- 
tions, which we have done. As we have mentioned the non-relativistic limit eq. 
(172) might become quantitatively incorrect in the vicinity of Ty "' rel yielding 
an overestimation of the density and, thus, an underestimation of T^ h . An- 
other reason might be that in reality g v is less than 10. For g v ~ 8 we would 
get higher values T^ Tcl ~ 140 MeV, T" h ~ 130 MeV. Nevertheless, within 
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Fig. 4. Logarithm of the ratio of the fermion density to the corresponding Boltzmann 
quantity for g v = 10 and g v = 7, and the effective transverse vector boson mass as 
functions of the temperature in units of T^t- The vertical dash-line indicates the 
area of the HBC. 



a reasonable variation of parameters we always obtain that the blurring of 
the nucleon continuum and the HBC of u and p mesons could occur already 
for T < m w . Concluding this discussion we again stress that we considered a 
model problem and disregarded many relevant couplings that might yield an 
effective repulsion. 



The typical behavior of the logarithm of the ratio of the fermion density to the 
corresponding Boltzmann quantity as function of the temperature is shown in 
Fig. 4 for T < T£ b . Moreover Fig. 4 demonstrates the temperature dependence 
of the effective transverse vector boson mass. We see a huge enhancement of 
ln(yOy/pj ]B oi) and a drastic decrease of m* bt /m h in the vicinity of T = T£ u ,T" h . 
At the critical point the first order phase transition to the HBC state may 
occur. The value m^/rrif, does not reach zero at this point. 
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6 Heavy fermions and less massive pseudo-scalar bosons (pseudo- 
vector coupling) 



Let us now consider pseudo-scalar boson - spin \ fermion system interacting 
via the pseudo-vector coupling 



Lint = -ig P .vm b ^§7 5 f#. (185) 

In realistic case bosons are 7r + ,7r~,7r° and fermions are n and p. The nNN 
coupling is g p , v ~ 1. 



6.1 Low temperature limit 



Again we start with the simplest $ given by the first diagram (2). We calculate 
the value 



S° = V (q)A° f (p + q)Vo(-q)A° b (q) 

= N p . v g 2 pv mf[2pqi - (m f + j>)q 2 + q 2 ^}, (186) 

where N p _ v = 3 is the isospin degeneracy factor (tt + , tt~ , 7r°). Within the non- 
relativistic approximation for fermions replacing po ~ to/, q 2 ~ to|, using 
qo <^ rrif, |g|<m/, and retaining only zero components of vectors (dropping 
linear terms in 7, which do not further contribute) we simplify (186) as follows 

S° ~ N p . v gl v mfm f [q 2 (l + 7o ) - g 2 (l - 7 o)] ■ (187) 

With the help of eqs (187), (22), (24) we calculate the fermion 3-momentum 
distribution 



«/,(±)(p) - 



{e p + m f >y ) (I + <y )N p . v g 2 v tp 



2e p 8vr 2 UuV to 6 

where Iq v (x) is determined by eq. (141). 
Finally, the fermion particle and antiparticle densities are 



«Boi,(±), (188) 



P.f,(+) PM-) 



PBol,(+) PBol,(-) 
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This result essentially differs from eq. (47) for scalar neutral bosons, since 
Iq v and Iq s are different. The difference with the result for the vector-isospin 
vector boson system for the case of a warm hadron liquid, where Iq v term is 
suppressed (see eq. (145)), is hidden in the different coupling constants and 
spin-isospin dependent coefficients. 

For T ~ m n /2, the regime of the warm hadron liquid is realized, mj ^> T ^> 
ml/rrif. We find p s p ym = p sym = p s p ym = p s ^ m ~ 1.16^^, that yields rather 
small correction due to a moderate value of the coupling constant, g w N — 1 <C 
g aN ~ 10. Nevertheless this correction is of the order of that for a and uo 
mesons. The reason is that the pion is lighter than o and uj. Thereby, Iq v is 
enhanced, see eq. (143). Moreover for T ~ m^/2 for o and uo still the regime 
of the slightly heated hadron liquid occurs rather than the regime of the warm 
hadron liquid, as for the pion. 

Contribution of the correlation diagrams is negligible in the low temperature 
limit, as in the scalar and vector boson cases. 



6.2 High temperature limit 



6.2.1 Analytic solution for fermion Green functions in STL approximation 

We continue operate with the single diagram of $ (the first diagram (2)). Now 
we calculate the value 



S = V (q)A f (p)V (-q)A° b (q)A f (p) 
= H.v m b 2 (-pV + 2(pg)#) (Gf) 2 



-3gl v m^ 2 (m 2 q 2 (G*) 2 + 2m f (q 2 $ - qp ■ i)G*G*) , (190) 

entering eq. (29). We used eq. (35), where we have put G 3 = 0. Dropping 
terms linear in 7 and qp and using (35), we may rewrite the Dyson equation 
(32) for the fermion sub-system, as the system of coupled equations 



p 2 Gf - mjGf = 1 + [(Jf- - JDp 2 + 2pgJT](Gf f 

+(jr- jrK(Gf) 2 , 

G* = Gf + 2J?- V G?G§, (191) 
where now 

J i V = J j2nY J 27 r b '(+) n6 >(+)te°) + Ab,(-)n b , H (qo)] , (192) 
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J 2 V = J J 1^ [Ab,(+)n b , {+) (q ) + A 6) (_ ) n 6i (_)(gb)] . (193) 

Setting Gi ~ G 2 in quadratic terms ((G R ) 2 ) we arrive at the Dyson equation 

Gf = + G} ,fl [(p§ + 3mJ) - (Jf" - jnf](G?) 2 

+jr(cf) 2 . (194) 

This equation already allows for the analytic solution. 
6.2.2 Intensity of multiple scattering 

Within the quasiparticle approximation for bosons and antibosons the spec- 
trum of boson excitations becomes 

iuf ±) (q,T) ~m\ + (3 p . v (T)q 2 + p p . v (T)q 4 + 

(195) 

cf. subsection 6.3. Due to the p-wave nature of the interaction (185) in eq. (195) 
there appear the bare boson mass and a significantly corrected q 2 term. Even 
one may have /3 P . V < in some temperature interval. For cases of sufficiently 
small j3 p , v and for f3 p , v < one should keep in (195) the higher order terms, at 
least the g 4 term. 

Using this spectrum we may evaluate the intensity of the multiple scattering 
Jf' v . For P pv > 0, assuming it is not as small, dropping the term [3 p , v q A in 
(195), with the help of eq. (192) we find 



JP . V _ G9l f qWq J dq 5(q 2 - m 2 -q 2 -Re^\q 2 )) 
. H.v 7 Q 4 d\q\ 1 



/ 



" 2n 2 m 2 J ( m 2 + p pv ^f 2 exp [ (m 2 + p pv pf 2 /T \ - 1 ' 
and in limiting cases 



7rV T 4 
jp.v Hp.v 

lOpptfm 2 



J*» = -J&L- 2 , for T»m b , (197) 
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J! 



p.v 



V 2 



03/2 3/2^/2 1/2 -exp (-m 6 /T) , for T < m 6 . 



;i98) 



Compare these expressions with the corresponding eqs (161), (162) for the 
vector boson case. 

Using (193) we obtain (for f3 p , v > 0) 



q 2 d\q\ (m 2 b + f3 p . v q 2 ) 



=»2\l/2 



2vr 2 m b 7 exp [{m 2 + (3 p . v q 2 ) l/2 /t] - 1 ' 



and in limiting cases 



(199) 



jp.v jp.v 

2 ~~ - ~ q3/2 o ' 



10/3 P °; z m 2 



for T > m 6 , 



(200) 



o 2 T 3/2 1/2 

2 = 23/2^3/2^3/2 ex P (~ m b/T) , for T « m 6 . 



(201) 



For — > intensities of the multiple scattering anomalously increase. How- 
ever in reality they do not tend to infinity, if one incorporated the term fi v . v q A 
in eq. (195). 

The maximum value of the density that can be achieved in this regime can be 
estimated as (84), now using eq. (200): 



ULcXX. 

Pf,(±) 



23/2 7r 2j3/4 



23/4^7/2^42 



(202) 



As we shall see in subsection 6.3, the quantity f3 p . v decreases with the increase 
of the temperature. For T = T^ h \ the value (3 P , V reaches zero and for T > T^ h \ 
it becomes negative. In refs [26,33] the density/temperature, when (3 p v crosses 
zero, was called the critical density/temperature for the appearance of the 
liquid phase of the pion condensation. There is no yet a long-ranged order for 
T^ h v > T > T^ h \. However there arise many virtual boson excitations carrying 
finite momentum \q c \. Actual values of particle momenta are near \q c \, and 
directions of the momenta are randomly distributed. For f3 p . v (T) < 0, the 
quantity 



ti 2 (q, T) = m 2 b + ReSf (q = 0, q, T), 



(203) 
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that can be called the effective boson gap, gets the minimum, cf. [14,33]. At low 
go, the boson spectrum has a non-quasiparticle nature. As we will demonstrate 
it below, see eqs (222), (224), the dispersion relation is then given by 



ifa(\q\)u ~ uj 2 +p (\q\ - \q c \) 2 , \q c \ ^ 0, u 2 , fa, (3 > 0, (204) 
where 

cu 2 = u 2 {q c ) =m 2 b + ReSf (q = 0, q = q c ) (205) 
The value 



-P P .v/(2P P . v ) (206) 

corresponds to the minimum of £; 2 (|g|), /3 — —2f3 p _ v . The form of the spectrum 
(204) coincides with that previously used in the description of the dense baryon 
matter in the pion condensation problem, cf. [14]. The boson spectral function 
has the form 



2, ( 207) 

Replacing (207) in eqs (192), (193) and using for simplicity the limit faT ^> uj 2 , 
P 1 = Pi(q c ), we calculate 



Ji v * o T%- , JV^Jr- (208) 

One has <C Jf' v , since J%' v has no singularity for uj q — > 0. We find an 
anomalous increase of the intensity of multiple scattering for u — > 0, i.e. the 
critical opalescence. In the given case, as for vector bosons, it is the precursor 
of the first order phase transition to the HBC However in difference with the 
vector boson case here the phase transition may occur to the crystalline or 
liquid crystalline state, see further discussion in subsection 6.3. 



6. 2. 3 Non-relativistic fermions 

In non-relativistic approximation for fermions we may put \p\ <C p and p ~ 
rrif. Then (194) simplifies as 

Gf = G° f ' R + G° f 'Hm 2 f J?- v (Gf) 2 . (209) 
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Its solution takes the form 



ReGf ~ ReGf ~ t^L^ ^ = g^m 2 , 



ImGf ~ ImG* ~ r_ >_ _ ; ( 2 10) 



2^ v -(p 2 -m 2 +i0) 



Ip.v 

for 2^f 2 v > (p 2 — m 2 ) 2 . The only difference with eq. (72) (the latter equation 
is valid for the scalar boson) is that in (210) enters j p , v instead of j s . There- 
fore, after the replacement J s — > Jf ' v expressions (73) - (83) hold also in the 
pseudo-scalar case. 

The temperature of the blurring of the fermion continuum evaluated within the 
non-relativistic approximation for fermions (Ty f' n ' rel ) follows from the relation 
Jf' v = m 2 /A (for /if = 0), from where using (197) (at the condition f3 p . v > 0) 
we estimate, cf. [5], 



6.2.4 Relativistic fermions 

Let us consider the contribution of the energy region p <C rrif. Then from 
(A. 19), (88), (194) we find the fermion- ant if ermion density (for one species of 
fermion) : 



T 2 I* 1 



8pu±)(po$T) = — J fd\p\ 



yjl2m)JZ v - AJf v p 2 - U v 2 v m) - (p 2 + m}) 2 
X 2m}[3Jr - J P 2 V - Jrp/m 2 ] ' ^ 

pL, = -m) - 2 jr + j?-y + ^jrm} - u$v f , 

cf. (89). We see that for — \ J% V > Jp}„ = (positive square root in 
(212) for p ~ 0) the fermion sub-system represents a rather dense packing of 
fermion- antif ermion pairs. The condition 



J v . v ^jr -\jT = jfi = m)l\2 (213) 
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determines the typical temperature T^ v f of the blurring of the fermion contin- 
uum. Assuming f3 PmV > 0, and using (213) and high temperature limit estimates 
(197) and (200), we obtain 



2 

rp.tnbl.f ^ m f 



* W^iU- (214) 



and 



bLf (S-^) 1 / 4 2V4 7r i/2^ (3 _^ ) i/4- ^ ^ 

For pions we find T^ v { ~ 215 MeV for (3 p v = 1. As a more realistic estimation, 
taking (3 P , V ~ 0.5 we obtain 1*£ ~ 132 MeV. 

For T > Ty' f, small energies po ~ T essentially contribute to the filling of the 
particle and antiparticle Fermi seas. For T in the vicinity of T^ v { , using that 
|Pmax| = for J p , v = Jf v l and expanding |p max | in < J p . w - J^ f < J™/, from 
(212) we obtain 



and 



(36 - 26/3 p ., + 4^jV 2 , 



(3-/? p .,) 7 / 2 t 2 [jr-(jn bl - fl2 



(216) 



~ T) ^ 3(7-2^)3/ M >~ • (217) 



For artificially large values J{' v ', Jf ' 1 " 3> m 2 from (212) we evaluate 



P/,(±)(PD^r)^ r ^. W) y , Plx^(3-/9 p .,)mj. (218) 
Using (197) we find 



10V 2 ^/> 2 m b (3 - p pM ) 
P/,( ± )(po S T) ^ ^_ . (219) 

If (3 P . V were independent of T, we would get a saturation of P/,(±) with in- 
crease of T. However j3 p _ v decreases with the temperature increase, see next 
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subsection. Maximum available density of fermions can be estimated equating 

(217) plus (79) (after the replacement J s — > Jf' v ) and, on the other hand, 

(218) . With the further increase of the temperature, (3 P _ V (T) diminishes and 
Pf,(±)(j>o ~ T) decreases. Eqs (217), (218) continue to hold for (3 p/u < 0, if one 
sets in these equations (3 P . V = (since J%' v is suppressed for f3 p . v < 0). For J{' v 
then one should use the estimate (208). 



6.3 Medium modification of boson excitations. HBC 



Now we are at the position to evaluate the boson self-energy. Taking the trace 
in (14) we obtain 



TrV (q)(p + m f )V (-q)(p + m f ) = -J^tf^ + mj) - 2p 2 Q q 2 ], (220) 

m b 

where we dropped small terms oc pq. Further for the sake of simplicity we 
will use the non-relativistic approximation for fermions and we put G± ~ G2, 
G3 = 0. Also, restricting ourselves by consideration of small g M , let us for 
simplicity omit a q dependence of the fermion Green function yielding the 
wave function renormalization terms. Then from (14) we find 



Re^(q)c-32N bar gl v mfq 2 m} J J^ReG«(p) 

^ ' 

x lmGK(p)[n fM (p ) + n/,(-)(po)]. (221) 

Nf, ar = 2 takes into account two type of fermions (neutrons and protons. Quite 
similar to that how we derived eq. (94) we obtain 



2N har gl v q 2 p f j 



m; 



1 - 



T 



p.v 



I2 1 



T 



■p.v 



(222) 



Here Pf j is, as above, the density of fermion-antifermion pairs for one fermion 
species. With this expression we are able to recover the quantity 



Pp.v 1 



2N bar g pAJ p }J 



(223) 
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entering the boson spectrum (195). This equation becomes invalid for rather 
small f3p mV , since we did not take into account q 4 terms, see (195). 

Instead of (109), that we had for the scalar boson, for the given case of the 
pseudo-scalar boson we obtain 




(224) 



Eq. (94 ) (after dropping there the wave function renormalization terms) and 
eq. (222), and also (109) and (224) transform to each other with the help of 
the replacements g 2 — > N^ ar g 2 v q 2 jm\ and J s — > Jf' v . From (224) for q Q <^ T 
we find ImSf ~ —Pi(q)qo, where the value 

&® = J,^WyW2 > fOT 4 V^T»^ ( 22 5) 

controls the low-energy part of the spectrum (204). In the boson self-energy 
for f3 p . v > we retained the term oc q 2 and dropped terms of the higher order 
oc q A . In order to find an appropriate spectrum for small f3 p _ v > and for 
j3 p , v < one needs to incorporate oc q 4 terms. 

The second order phase transition to the p-wave HBC state could occur for 
f3 p .v < 0, if the effective boson gap ujq reached zero at some finite momentum 
q c . However the eq. (205) has no solution for ojq = at least in the region of 
the validity of the non-relativistic approximation for nucleons. This is due to 
the fact that the intensity of the multiple scattering entering equations for the 
density pjj and ujq is Jf' v oc l/a>o, see eqs (205), (208) and (222). The second 
order phase transition is, thus, excluded. The phase transition to the p-wave 
HBC state is then the first order transition. The latter possibility, however, 
needs a more detailed analysis (including the comparison of thermodynamic 
potentials of two phases). It is worthwhile to mention that the pion condensa- 
tion in dense nuclear medium also arises by the first order phase transition at 
finite temperature, namely as the consequence of thermal fluctuations [14,28]. 
The latter are described by the same tadpole diagram, as the intensity of the 
multiple scattering discussed here. 

The typical behavior of the logarithm of the ratio of the fermion density to 
the corresponding Boltzmann quantity as function of the temperature is shown 
in Fig. 5 for T < T^ h v . Here the effective boson mass does not significantly 
deviate from the bare mass but the p-wave boson-fermion interaction changes 
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Fig. 5. Logarithm of the ratio of the fermion density to the corresponding Boltzmann 
quantity for g p _ v = 1 as function of the temperature in units of T^i. The vertical 
dash- line indicates the area of the HBC. 

substantially q 2 terms in the spectrum. As in scalar and vector boson cases, 
we see a huge enhancement of ln(pj/pj B oi) m the vicinity of T = T£ u . At the 
critical point the first order phase transition to the HBC state may occur. 



6.4 Correlation effects 



The first diagram of $, that we have studied, sums up all perturbative dia- 
grams for the fermion Green function without crossing of the boson lines, cf. 
eq. (8). For the fermion-boson coupling under discussion, within the STL ap- 
proximation, each diagram with the one crossed boson line (see eq. (9)) brings 
the suppression factor v = v T v a compared to the corresponding diagram with 
the same number of boson lines but without their crossing, cf. [5,26]. Indeed, 
in the non-relativistic approximation for nucleons, v T = 1/3 is due to the non- 
commutation of the isospin r matrices and in non-relativistic approximation 
for nucleons v a = 1/3 is due to the non-commutation of the spin a matrices. 
Each product t% Tj Tj tj yields factor 9, whereas each product TiTjTiTj yields fac- 
tor 3. The same statement is valid for a matrices. Concluding, we may retain 
only the first diagram (2) in $ with an appropriate accuracy. 
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7 Towards the description of the state of hadron porridge 



We discussed the behavior of model systems with a small fermion chemical 
potential, assuming for the sake of simplicity only one type of the interaction 
in each case, like the case of spin | fermions (e.g., N) coupled with the scalar 
boson (<r), or with the vector boson (u or p), or pseudo-scalar boson (n). 
This is, of course, a gross oversimplification. In reality nucleons couple with 
all mesons: a, n, u, p, K, etc. Also high lying baryon resonances, like A- 
isobar, A/"* (1440), hyperons, etc., interact with the nucleon and with each 
other. Moreover, mesons interact with each other. We may call such a system 
state the hadron porridge, bearing in mind that already at sufficiently low 
temperatures the quasiparticle approximation fails to describe baryons and that 
for T > m n the continuum is blurred for all relevant hadrons, as we have shown 
it above on different concrete examples. 



7. 1 Hadron porridge in the low temperature limit 



In the low temperature limit contributions of different mesons are summed 
up yielding resulting baryon resonance distribution and the total density of 
baryon resonances. The density of the spin | baryon resonance of given species 
is easily recovered with the help of eq. (119): 



fm B *T\ 3/2 I" mg. =F p, B * 
Pb*,{±) = N B . [~^- J exp y 

Ar _/m B »T\ 3/2 
+N B « ]T ' 



2tt 



exp 



m B * =F Hb* 



(226) 



T 



" 2 ( m B*\ o 2 tO ( mB \ 

X 9aB*B*-L0s I - I + ^9ttB*B*-^0v [ ' 

+9„B*B* \1 0V [—) + 2— J 0v [ — 

^ i 

+3o 2 (jo 

+69 ^{ Iov {m p ) + 2m Bt Iov 



4tt 2 




The summation is performed over all possible s B * = 1/2 states including 
the given baryon state {B*). Dots mean contributions of other meson ex- 
changes and transitions between baryon states with different spins in initial 
and intermediate states in the diagram (6), e.g. (N — K — A; A — K — N), 
(N — it — A(|, |); A(|, |) — 7r — N), etc. Correlation effects are unimportant 
in the low temperature limit. 

Numerically, already the contribution given by terms, which are explicitly 
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presented in eq. (226), yields Rn = Pn/Pb \ — 1-5 for T ~ m w /2 and for 
values of coupling constants, which we have used above; pg ol is the nucleon 
density calculated with the help of the Boltzmann distribution. With taking 
into account of other contributions one may expect a further enhancement of 
the ratio. 

We would like to pay attention to the enhancement of the population of baryon 
resonances, e.g., of the strange baryons, since contributions like (A—tt—N; N— 
7r — A), (A — K — N; N — K — A), etc, have no any additional exponential 
suppression compared to (N — it — N; N — n — N) terms, see eq. (120) and 
the corresponding discussion. Thus, the ratio of the densities of the given 
resonance and the nucleon, Rb*/n — Pb*/pn, can be ~ 1 for the low lying 
baryon resonances, like A, A, etc, for which B* -meson-N coupling constants 
are not suppressed. 



7.2 Hadron porridge in the high temperature limit 

Description of the hadron porridge in the high temperature limit, T > m n , is 
much more involved. Even in the STL approximation one gets coupled system 
of Dyson equations for different baryon resonance Green functions. Intensities 
of the multiple scattering entering the Dyson equation for the given resonance 
describe the coupling of this resonance to other resonances via the radiation 
and the absorption of different type of mesons. Correlation effects contribute 
to each vertex. The main correlation effect comes from scalar meson and vector 
meson terms. Moreover, boson effective masses and the \q\ dependent terms of 
the boson spectrum are changed. In addition to all has been said the problem 
should be solved self-consistently. 

Assume for a moment that correlation effects are suppressed. Then one deals 
with the sum of one-diagram $-terms, see first diagram (2), where the bo- 
son line corresponds to any relevant meson and the fermion line relates to 
one of baryon resonances. Summation is over mesons and baryon resonances, 
including nucleons. Assume also that baryon resonance states decouple, i.e., 
coupling constants g^B*,B* = 0, for % ^ j, /j, is one of mesons. Then, above re- 
sults obtained for the system with one type of baryons and one type of mesons 
are easily generalized to include all above mentioned types of the interaction. 
The contribution of the given baryon resonance B* in non-relativistic approx- 
imation for baryon resonances is obtained with the help of the replacement of 
the corresponding J to the sum of relevant J's, e.g. 

J s -> J*,b* + Ji' B * + J P i' B * + JT B * + - (227) 
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Inclusion of cross-channel couplings (#»£*, s* ^ 0) and correlations make the 
problem much more involved. Correlations may yield a significant suppression 
of the resulting intensity of the multiple scattering compared to that given by 
(227). Moreover, going beyond the non-relativistic ansatz for baryons, that 
is necessary to do for T > T b i.f, additionally complicates the consideration. 
We postpone a more detailed discussion of properties of the hadron porridge 
(including baryon resonances with a cross-channel coupling) to the subsequent 
paper. 



8 Conclusion 

The description of the hadron system with a small baryon chemical potential, 
|/-*bar| <C T, at the finite temperature is essentially different from the descrip- 
tion of the dense baryonic system |/^bar — ™n\ 3> T, cf. [14]. The former regime 
might be relevant for heavy ion collisions at RHIC and LHC, whereas the latter 
one, for SIS energies. 

We treated the fermion-boson problem self-consistently considering first one 
heavy fermion and one lighter boson species. Then we estimated a role of 
correlation effects. 

• We found that already at low temperatures, T <C T hU ~ m n ~ 140 MeV, 
the nucleon 3-momentum distributions are enhanced compared to ordi- 
nary Boltzmann distributions. Quasiparticle approximation fails already in 
regime of a warm hadron liquid (ml/rrif < T <C minjmfe, Tbi.f}), if this limit 
is indeed realized for the given particle species and if the coupling constant 
is g > 1. 

• For a higher temperature (typical value is T ~ T b i.f ~ m^) the hadron 
continuum is totally blurred. Effective masses of all hadrons are decreased 
and baryon widths are dramatically increased due to baryon rescatterings 
on bosons. Baryons play a role of a glue for mesons. The hadron liquid 
comes to the state of a hadron porridge. 

• Scalar mesons, like a, may undergo second order phase transition to the 
s-wave hot Bose condensate state for T > T c s b ~ T b i.f ~ m n . 

• Vector mesons, like u, p, may undergo first order phase transitions to s- 
wave hot Bose condensate states for T > T" h ~ Tbi.f — Tn n . in the vicinity 
of the transition point the intensity of the multiple scattering is significantly 
increased leading to a critical opalescence. 

• Pions may form a liquid hot Bose condensate for T > ~ m n and with 
a subsequent increase of the temperature may undergo first order phase 
transition to the p-wave hot Bose condensate state. As in the vector case the 
intensity of the multiple scattering is significantly increased in the vicinity 
of the critical point. 
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• For T > Tbi.f fermions fill a rather dense packing state. 

• The number of baryon resonances, like A, should be not much less than the 
number of nucleons. In particular cases, like for A isobar, one even may 
expect a higher abundance than for nucleons. The number of kaons, which 
interact with nucleons and hyperons, should be also large. One may expect 
the same effects for them, as for pions and cr's, due to a strong s- and p-wave 
KNA interactions. 

When scalar meson mass reaches zero, the alternative state to the hot Bosc 
condensate one could be the chiral phase state. To conclude, which state is 
preferable, one needs to study the chiral model, e.g., the linear sigma model, 
and one should include quark effects. 

For the low temperature regime, as a guide, we used the Urbach rule. As a rel- 
evant approximation to simplify the self-consistent analysis of the hot fermion- 
boson system (T > m n ), we exploited the soft thermal loop approximation. The 
latter approximation allowed us to solve the problem. The self-consistency of 
the approach is manifested in the fact that fermion Green functions completely 
loose their former pole shape. It is the consequence of multiple collisions of the 
fermion with bosons. On the other hand, bosons drop their masses interacting 
with modified fermions, that further stimulates the enhancement of the fermion 
width. Please notice that with only one $ diagram we estimated Tbi.f < m 7r . 
Attractive interactions of the nucleon with other meson species, being taken 
into account, could still decrease the value T bLf . However we have shown that 
correlation effects in the scalar boson channel increase the value T blf . In the 
regime that we called "the hadron porridge" interactions between many reso- 
nance and meson species come into play. Scalar bosons contribute to vertices 
with participation of other mesons effectively suppressing vertices. We consid- 
ered only some of interactions. Therefore it is hard to realistically estimate the 
actual value of T bbf . Our rough estimates show that T bLf is somewhere in the 
vicinity of the pion mass (140 MeV). However further more detailed analysis 
is necessary. 

We described the system in terms of only hadron degrees of freedom and we 
found a state of a hadron porridge for T > m n . In realistic situation quark- 
gluon degrees of freedom are also excited for T > m n . These effects may yield 
an additional diminishing of effective hadron masses due to the decrease of the 
quark condensate density with the increase of the temperature. The role of 
quark- gluon fluctuations rises with increase of the temperature. More likely, the 
system comes to a strongly correlated hadron- quark- gluon state or, in another 
words, a boiled hadron- quark- gluon porridge, (HQGP). Deconfinement, in a 
standard meaning, as the pure quark-gluon state, is probably delayed up to a 
higher temperature. 

Our aim was to discuss qualitative effects. Thereby, we focused our study on 
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simplified models. In order to quantify the consideration and to do an ap- 
propriate fit of experimental results, cf. [1,2,3], one needs to include A and 
other nucleon resonances and hyperons. In the meson sector one has to in- 
corporate kaons. Meson-meson interaction, short-range correlations and other 
effects may also contribute essentially, as we have argued. Thermodynamical 
characteristics, like the pressure, energy, entropy still should be computed. 
Comparison with lattice results, cf. [12,34], is then required. Moreover, one 
should supplement the consideration of in-medium properties of the system 
by a relevant model of the freeze out stage to match in-medium particle dis- 
tributions to those distributions one observes at infinity. For sure the realistic 
consideration of strongly interacting heated nuclear system is a very compli- 
cated problem that needs a further study. We hope to return to these questions 
in forthcoming publications. 
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APPENDICES 
A Formalism 

A.l Green functions and self-energies 

In matrix {— , +} notations there exist exact relations between two-point func- 
tions F (Green functions and self-energies), cf. [4], 

pR = p— _ p-+ = p+- _ p++^ 

p A = p— - p+- = p~+ - (A.l) 

F R and F A are retarded and advanced two-point functions (operators in spin 
space). The co-contour functions are Fjk = <Jji<Jk m F lm , a* = 5*, = are 
standard 03 Pauli matrices, space indices j, k,l,m — 1, 2, 3, 6* — 1 for j = k 
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and 5j = for j ^ k. 2 

The free fermion retarded Green function is 



= p2_ A ly +i0 - G / R (9)A?W. (A.2) 
and the free boson retarded Green function is 

A° and A° are corresponding spin structure operators of fermions and bosons, 
p, q are 4- momenta. For the spin 1/2 fermions one has A° = $ + rrif and for 
the scalar neutral boson one has A° = 1, $ = 7^q^, 7 M are Dirac matrices, 
/i = 0,l,2,3. 

The spin structure of the free massive vector-boson propagator is given by 

A° b = (g^-q^/m 2 b ). (A.4) 

Using (1) it is convenient to introduce one-particle spectral and width func- 
tions (operators), 



Aj = -2ImGf (q) = -21m ^ 1 ^ , f < = -2ImEf . (A.5) 

M; + Ty2 

The quantity 



Mj = (G^)" 1 - ReEf (A.6) 

demonstrates the deviation from the mass shell, M = on the mass shell in 
the matter, the relation that introduces quasiparticles. 

At equilibrium all non-equilibrium Green functions, G , G ++ , G h , and 
G + ~, are expressed through the retarded Green function, cf. [4], e.g., 



iG~+ = T n t (q )A t (q), iG+~ = (1 T n t (q ))Mq). (A.7) 



2 We use notations of ref. [4], where Green functions G y (here i,j = {+,—}) are 
the same, as in refs. [16,21], and S ^ £ H differ by the sign from quantities used 
in refs [16,21]. Self-energies £y coincide with S 4 - 5 of [16,21]. 



76 



The upper sign corresponds to fermions, the lower, to bosons, 



exp[(g - lk)/T\ ± 1 

are fermion/boson occupations, fa are chemical potentials. Relations for 
and are similar to eqs (A. 7): 

= = Tni(q )Ti(q), 

= = (1 T n, t (q ))f t (q). (A.9) 



A. 2 P article- antiparticle separation 



To deal with Green functions of physical particles and antiparticles we intro- 
duce the decomposition of Green functions into particle (+) and antiparticle 
(— ) parts 

Gf{q) = Q{qo)G? M {q) + 6(-g )G^ ) (-g), 
GT + (q) = &(qo)G~; +) (q) + e(-gb)G+- ,(-?) , (A.10) 



= e(g )^,(+)(g) - e(-g )A,(-)(-g), 



(A.11) 



where O(rr) = 1 for rr > and =0 for rr < 0. These equations allow to 

reduce the integration in q in self-energy diagrams to positive energies. 

Spectral functions satisfy sum rules: 



oo 

Nf J 



To 



Af,(+)(qo,q) +A fiH (qo,-qj\ ^ = 2, 



(A.12) 



^Tr 



qo 



A,(+)(qo, q) + A bi (-)(qo, -q) 



o 



2tt 



1. 



(A.13) 



The trace is taken over spin degrees of freedom. Nf = (2s f + 1) is the fermion 
degeneracy factor and N b = (2s 6 +l) is the boson degeneracy factor. For vector 
bosons sum rules are fulfilled separately for transversal (N b = 3, A bt (±) — > 
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-<4&(±)) and longitudinal (N b = 1, ^4&,(±) — > ^>(±)) degrees of freedom, see eqs 
(133), (134). 

At the thermal equilibrium we may introduce occupations of physical particles 
and antiparticles as 



"»,(+) (?o) = ni(q )Q(qo), n»,(-)(go) = -M~?o) =F l]0(5o), ( A - 14 ) 



/ x ©(5b) 
^,(+)(?o)- 



exp[(g -^)/T]±l' 

= ex P [(* + ft )/r]±r (A ' 15) 

In-medium 3-momentum distributions of particles and antiparticles are ob- 
tained with the help of (A. 15) from the definition of the zero component of 
the Noether 4-current, cf. [4], 



M _ ? ffiA(+)(g) d Qo 
n ^ +)[q) J exp[(9b-^)/T]±127r' 



oo 



(-)(?) = / f Al ' ( i q \ ^, (A.16) 
i,( )W J exp[(g + /i,)/T] ± 1 2vr' 1 J 

£° = 70 for fermions and vf = 2q for bosons. The chemical potential of the 
antiparticle differs by the sign from the chemical potential of the particle. For 
the non-strange bosons the boson number is not conserved and the chemical 
potential fib = 0. 

For a large size system we may neglect particle radiation from the surface 
and the baryon chemical potential is determined using the baryon number 
conservation 



Pbar = P/,(+) - Pf,(-)- (A. 17) 

The total Noether 4-current presented as the sum of partial contributions 
is conserved quantity [4]. It is expressed through the standard Ai spectral 
functions. As it has been shown in [35], in thermal equilibrium the Noether 
current coincides with the conserved effective current (that can be derived 
from the Kadanoff-Baym kinetic equation), provided $ derivable approxima- 
tions are used. Refs [4] have demonstrated thermodynamical consistency of 
the approach. Thus the zero component of the Noether current coincides with 
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the corresponding quantity obtained with the help of the standard thermody- 
namical relation, e.g. for the baryon density 



Pba 



bar 



(A.18) 



Using the latter expression ref. [36] studied the nucleon - delta isobar - pion 
system in the virial limit of Dashen, Ma and Bernstein [37]. 

Particle and antiparticle Noether densities are 



Pi,(±) = Tr / "i,(±)(9) 



d 3 q 

W) 



3 ' 



(A.19) 



In the paper body we focus on the discussion of only slightly particle-antiparticle 
asymmetric systems and, as the limiting case, the particle-antiparticle sym- 
metric system (heated hadron vacuum), bearing in mind possible applications 
to heavy ion collisions at RHIC and LHC energies, respectively. We assume 
Pbar "C Pf,{±) and, thereby, \pi\/T <C 1. Then one may simplify above equations 
as 



n 



<:,(+) (<?o) 



n 



fob) 



l + ^(l=Fnr(9b)) 



(A.20) 



"»,(-) (9o) = rii ym (%) 



Pi 



(A.21) 



where the fermion/boson quantities are supplied by index "sym" for the 
particle-antiparticle symmetric system. For p bar > p/,(_) particle-hole effects 
become important. 



B Non-relativistic heavy fermions 



In the paper body we have shown that in the low temperature limit heavy 
fermions can be treated as non-relativistic particles, since both the energy 
deviation from the bare fermion mass and the fermion momentum are much 
smaller than to/. In the high temperature limit the non-relativistic approxima- 
tion for fermions also holds in some temperature interval, namely, if conditions 
(70), (165) are fulfilled, the former for the case of the scalar or pseudo-scalar 
boson exchange, and the latter for the vector boson exchange. The relativistic 
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study performed in the paper body justifies the validity of the non-relativistic 
approach for given cases. In order to come to non-relativistic expressions from 
relativistic ones we need to expand G\ ~ G<i near the fermion mass-shell, 
p 2 = mj, and to introduce the fermion non-relativistic self-energy (38) aver- 
aged over spin degrees of freedom. 

Suppose, conditions for the non-relativistic approximation for fermions are ful- 
filled. Then the consideration is simplified. The free fermion particle/antiparticle 
Green function obeys the equation 



[G°/ R ( P0 ,p)r cl = r~ ~ , (B.l) 

instead of eq. (A. 2) in relativistic case. The spectral function (instead of (A. 5), 
(A. 6)) is given by 



pn.rel 

/jn.rel _ / (U O) 

f ~ [Mf™ l } 2 + [r n - rcl ] 2 /4' 1 ' ; 

with Mf- rcl = Po~p 2 / (2m/) — [ReSj] n - rel . The sum rule is changed accordingly. 
Fermion particle and antiparticle distributions and densities are 



nV el (q)= / n 1 \ , n — , (B.3) 

/ W / exp[(g T/i/)/T] + 127r' 1 ' 



cf. eqs (A. 16) and (A. 19). Further to shorten notations we suppress index 
"non — rel". 

Thermodynamic quantities are easily found. The Noether energy-momentum 
tensor [4] 



is exactly conserved within the $ derivable scheme. Here V is the system 
volume, i = {/, b} indicates that sum is over all non-relativistic fermions and 
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relativistic bosons. Extra factor (j^J occurs only for neutral bosons. The 

velocity v? = dG^/dp^. The potential energy £ pot , which a probe particle 
with Wigner density G~ + would have due to the interaction with all other 
particles in the system, is 



£ pot = Tr 



2)n.b.?Z 



d A pi 



R 



(B.6) 



and the interaction energy £ mt is simply related to £ pot , 



gint _£P ot 

a 



(B.7) 



in case all the interaction vertices of the theory have the same number of lines 
attached to them. For the two-fermion-one-boson interaction, see (2), one has 
a — 3, that results in £ mt = |£ pot . The additional term £^ CT ^ appears in 
eq. (B.5) only in the case of derivative coupling, cf. [35]. The energy and the 
pressure at the quasi-equilibrium are given by 



E 



e 00 , p = -n/v = - (e 11 + e 22 + e 33 ) , 



(B.8) 



n = VTr £ [ -^n t (p 0t ) {-21m In \-Gf ( Poi + i0,#) 



(B.9) 



$ is the functional expressed in terms of non-equilibrium contour Green func- 
tions, cf. [4] . The expression for the entropy is as follows 



S ° = S/V = -Tr]T J 



d% dni(p 0i ) 



(2tt) 4 dT 



x 



-21m In [-G?(poi + iO,#) 



- ReGf r, - AReE- 



(B.10) 



This way one may recover all necessary thermodynamic characteristics of the 
system. 
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